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Abstract—The effect of small additions of sodium oley|-p-anisidinesulphonate on the rate of 
extraction of iodine from aqueous solutions with 4-2 to 22:3 mm? falling droplets of carbon 
tetrachloride has been determined. The droplets were formed in a section of the column con- 
taining water and surface-active agent only, thus eliminating extraction during drop formation 
and securing a stationary adsorption state before the droplets entered the extraction column 


proper. 

A maximum reduction in the overall mass transfer coefficient of 67 to 68% was observed for 
all droplet-sizes, with an addition of about 6 x 10 g surface-active agent per 100 ml aqueous 
solution, whereas much larger additions were practically without further effect. 

The results have been analysed in terms of an interfacial resistance coefficient R, sec/cm, 
which was found to be independent of the droplet size, and to obey the equation Cy = 5-4 x 10% 

Ry 


where Cg is the aqueous concentration of surface-active agent in g/100 ml, and R, 


is calculated on a solvent basis. The equation is derived on the assumptions that the interfacial 
resistance is proportional to the interfacial concentration of surface-active agent, and that the 
adsorbed film is of the “ vapour expanded" type. An adsorption isotherm for a “ vapour 
expanded " film is given. 

It is concluded that the results cannot easily be explained in terms of interfacial turbulence 
or transfer of turbulence, but that there is some kind of interaction between the adsorbed film 
and the diffusing solute, as postulated by Hutcuinson [2] and Garner and Hace (3). 

The rate of fall of the droplets was also found to decrease sharply on the addition of small 
amounts of surface-active agent, the reduction amounting to 7 to 10% for all droplet sizes at 
about 2 to 8 x 105g per 100 ml. Much larger additions resulted in a more gradual decrease. 


Résumé—Les auteurs ont étudié l'influence de petites additions d’oleyl-p anisidine sulfonate 
de sodium sur la vitesse d’extraction de l'iode & partir de solutions aqueuses dans lesquelles 
tombent des gouttelettes de CCl, de 4,2 & 22,8 mm 3 formées dans une section de la colonne 
contenant uniquement de l'eau et l'agent tensio-actif. On élimine ainsi toute extraction pendant 
leur formation et l'on assure dans la colonne d'extraction proprement dite état stationnaire 
avant que les gouttelettes n'y entrent. 

Dans le cas d'une addition d’environ 6-105 g. d’agent tensio-actif pour 100 moles de solution 
aqueuse, la réduction maximum du coefficient global de transfert de masse est de 67 a 68% 
quelles que soient les dimensions des gouttelettes ; des additions plus importantes n'ont pratique- 
ment pas d’effet supplémentaire. 

Le coefficient de résistance 4 l'interface R;/sec/em est indépendant de la dimension des 


gouttes et obéit a I'équation 
Cs = 5.4x 
(1 


"© A summary of this paper has appeared elsewhere [1]. 
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ot Cg est la concentration aqueuse de l'agent tensio-actif en g/100 ml et R, est calculé sur la 


Les auteurs supposent qu’é l'interface la résistance est proportionnelle & la concentration 
en agent tensio-actif et que le film adsorbé est du type “ vapeur détendue.” Ils donnent pour 


ce cas, un isotherme d’adsorption. 


Tis concluent que les résultats s‘expliquent difficilement par une turbulence a l’interface 
ou un transfert de turbulence mais qu'il y a plusieurs types d’interaction entre le film adsorbé 
et le soluté qui diffuse comme |l'ont supposé Hutchinson, Garner et Hale. 

La vitesse de chute des gouttes décroit nettement par addition de petites quantités d'agent 
tensio-actif, de 7 A 10%, quelles que soient leurs dimensions, pour environ 243 x 10-5 g/100 ml 


INTRODUCTION 


Increasing attention is being devoted to inter- 
facial resistance to mass transfer processes, and 
the evidence collected indicates that in liquid- 
liquid extraction the Lewis and Whitman two-film 
concept (which forms the basis of conventional 
design methods) is not universally applicable. In 
addition to molecular and eddy diffusion, there 
is apparently in some cases a further resistance 
closely associated with the interface itself. This 
contention has been supported mainly by obser- 
vations that additions of surface-active agents in 
some cases materially reduce the rate of extrac- 
tion. Some investigators, however, have failed 
to find this effect, and Hutcuinson [2] has 
reported that adsorbed films which retarded some 
species of diffusing molecules did not retard 
others. The most important papers on the effect 
of surface-active agents on the rate of extraction 
have been reviewed recently by Garner and 
Hae [8], who have themselves provided strong 
evidence for the existence of interfacial resistance. 
Murpocn and Pratt [4] have shown that the 
rate of transfer of uranyl nitrate from aqueous 
solutions to methyl-i-buty] ketone and dibutoxy- 
diethyl ether is partly determined by an inter- 
facial resistance caused by a slow chemical 
reaction with the solvent molecules at the 
interface. The recent work of Lewis [14] with a 
new type of transfer cell will be discussed later. 
SrvuKE [6] has reported that some surface-active 
substances greatly reduce the rate of rise of oxygen 
bubbles through water. With increasing equivalent 
diameter of the bubbles, the rate of rise normally 
goes through a maximum at about 1-7 mm at 18°C, 
but with an addition of 5-2 x 10 g caproic acid 
per 100 ml this maximum disappears, and the rate 


dagent tensio-actif. Des additions plus importantes entrainent une diminution plus lente. 


of rise is reduced to less than half for this bubble 
diameter. Formic acid, on the other hand, had no 
such effect, even in concentrations which caused a 
twenty times greater reduction in surface tension. 

In a previous paper Hotm and TerseseEn [1, 5] 
have shown that additions of 8 x 10° g per 100 
ml of sodium oleyl-p-anisidinesulphonate reduced 
by about 40% the rate of transfer of iodine from 
an aqueous solution to carbon tetrachloride in a 
stirred extractor. These results pointed clearly to 
a surface resistance, but they were not susceptible 
to quantitative treatment, because the observed 
changes in the rates of extraction resulted from 
two opposing effects. The surface-active agent 
caused an increase in the interfacial area as well 
as a retardation of the rate of mass transfer per 
unit area, the latter effect predominating at con- 
centrations below 8 x 10° g per 100 ml, and the 
former above. 

In the present investigation, falling drops of 
carbon tetrachloride of known size were used. 
With the extremely low concentrations of surface- 
active agent, it is not safe to assume that a 
stationary adsorption state is rapidly established, 
and methods of correction for extraction during 
drop formation would be particularly uncertain. 
These difficulties have been overcome by forming 
the drops in a separate section of the column 
equipped with a rod type stirrer, and containing 
water and surface-active agent only. The drops 
were allowed to descend in a spiral path through 
this liquid under the action of the agitator before 
reaching the extraction section. It was also 
considered essential to prevent the liquids from 
coming into contact with rubber or similar 
materials. The rate of fall of the droplets in the 
extraction section was also noted in order to 
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ascertain any correlation between the effect of 
the surface-active agent on the rate of mass 
transfer and on the rate of fall. 

In aqueous solutions of iodine, the surface- 
active agent was found to add two atoms of 
iodine, presumably at the double bond. It is 
therefore the effect of this iodinated compound 
which has been observed in these investigations. 


APPARATUS AND PROCEDURE 
The apparatus is shown diagrammatically in 
Fig. 1. The upper section, called the adsorption 
column, was provided with a rotating glass rod, 


Overflow siphon 


burette 
25 mi. 


Fig. 1. 


extraction apparatus. 
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and contained an aqueous solution of surface-active 
agent of the same concentration as used in the 
extraction column underneath. The bottom end 
of the adsorption column was gradually narrowed 
down to an internal diameter of 7 mm, and ended 
in a Pyrex 10/30 ground-glass joint fitted with a 
hollow plug. This plug could be pushed out by 
carefully lowering the glass rod, and would then 
float on the liquid in the extraction column. When 
the glass rod again was raised, an equivalent 
volume of water containing surface-active agent 
was added from a burette at the top, to prevent 
the iodine solution in the extraction column from 
penetrating into the adsorption column. The 
extraction column, which had an internal dia- 
meter of 50 mm, was available in three lengths, 
giving alternative extraction paths of 1430-13810, 
940-830 and 580-440 mm, the ranges being due 
to the varying liquid level in the 25 ml burettes 
attached to the bottom of these columns. Both 
columns were provided with jackets for circulating 
water which was thermostatically controlled at 
25-0°C, Accurate temperature control was 
essential for preventing thermal convection 
currents, which would cause mixing of the liquids 
in the two columns. One hour was allowed for 
establishing thermal equilibrium before the plug 
separating the two liquids was pushed out. 

All rubber connections were so arranged that 
a cushion of air prevented the liquids from coming 
into contact with the rubber. The stop-cock on 
the burette was lubricated with glycerol and 
bentonite, whereas the one in the carbon tetra- 
chloride feed line was not lubricated at all. 

The arrangement for producing the drops can 
also been seen in Fig. 1. Carbon tetrachloride 
saturated with water was run from a constant 
head flask through a capillary which dipped into 
the liquid in the adsorption column. By adjusting 
the effective difference in level, and by using 
capillaries from 0-1 to 8 mm in diameter, the rate 
of formation of the drops was varied between 1 
and 8 per sec, and the size of the drops betwecn 
4-2 and 22-3mm*. Some difficulties were, how- 


ever, experienced in maintaining a constant drop 
size from one experiment to another. 

During un experiment the drops were finally 
collected in the burette at the bottom of the 
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extraction column, and the liquid withdrawn for 
analysis. The aqueous phase remained stationary 
in the column during each experiment, and was 
sampled before and after. Only a very small 
fraction of the solute was extracted. 

The drop volume was determined by counting 
the number of drops required to give a volume 
of 1 to 5ml in the burette. The rate of fall of 
the drops was obtained by measuring with a 
stop-watch the time for a fall of 1,000 mm, and 
the average of 10 observations taken. With 
varying rates of drop formation the distance 
between the drops varied from 70 to 500 mm, but 
this did not noticeably affect the rate of fall. 


MATERIALS 


The aqueous phase was not saturated with solvent 
before the experiments, as this was not found to 
influence the results. The solvent was, however, 
saturated with water. 

The iodine used was “ B.D.H. Laboratory 
Reagent, Resublimed.”’ For each experiment a 
small amount was ground in a mortar, and 
dissolved in water under vigorous stirring to give 
4 concentration between 180 and 210 mg iodine 
per litre, and finally filtered. 

Carbon tetrachloride Merck Puriss” was 
distilled in a packed column of 5-10 theoretical 
stages at a reflux ratio of about 10:1. The 
fraction boiling between 76.5 and 76-8°C was 
collected. Used solvent was regenerated by 
treatment with sodium thiosulphate followed by 
simple distillation. 

Sodium oley]-p-anisidinesulphonate was purified 
and dissolved in water, together with sodium 
sulphate and sodium chloride, as described 
previously [5]. 


Meruovs or ANALYSIS 


Solutions of iodine in water and carbon tetra- 
chloride were analysed photometrically using a 
Beckman DU Quartz Spectrophotometer, the 
readings being taken at the absorption maxima, 
which were found to be located at 4,600 and 
5,150 A respectively. The calibration curves 
giving optical density as a function of concen- 
tration were straight lines through the origin for 
both solvents. Samples were kept in glass- 


stoppered bottles, and analysed as quickly «s 
possible to avoid loss by evaporation. The 
formation of a compound between the iodine 
and the surface-active agent previously mentioned 
did not interfere with the analysis, due to the 
extremely low concentrations of surface-active 
agent used. 


EFrFrects 


Although the technique of forming the drops in 
solute-free aqueous phase eliminated the need to 
correct for extraction during the formation of the 
drops, the need to correct for extraction at the 
interface in the burette still remained. To determ- 
ine the magnitude of this correction, it was not 
considered sufficient to observe the uptake of 
iodine by a stagnant volume of solvent in the 
burette, because the turbulence caused by the 
drops striking the interface would entail a higher 
rate of transfer during actual opertaing condit- 
ions. As the nearest approach to determining 
true correction factors, a series of experiments 
was carried out with drops of equal size but with 
differing frequencies of formation. By comparing 
the observed mass transfers with the values extra- 
polated to an infinitely high frequency of format- 
ion, it was found that an average of 2-4 mg of 
iodine per minute were transferred at the inter- 
face. With the longest column the correction 
amounted to less than 1% with drops of 22 mm" 
volume formed at a rate of 3 per second. Most 
of the experiments reported in this paper were 
carried out under conditions where this correction 
was insignificant. 


RESULTS 


The effect of additions of surface-active agent on 
the rate of extraction was determined for three 
different droplet sizes, using the largest column 
where the droplets travelled an average distance 
of 1837cm. The experimental results have been 
summarized in Table 1. During a run, only 
between one and two per cent of the iodine in 
the aqueous phase was extracted, and it was 
sufficiently accurate to use an average concentr- 
ation for the calculation of the approach to 

Ce x 100 


equilibrium, p = a Prarce and Ever- 
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Results with 187 cm average extraction length. 


Table 1. 
CRO Cr 


mg/ml mg/ml 


Cg 


_ 
“mCra 


f 


Kew. 10 105 
Di cm/sec Dimen- | Ry, 10% 
sionless | sec /om 


Droplet volume : 22-8 mm® + 0-1, ar tae 


175mm + 0-005 


0 0-180 0-177 0-408 2-54 0-992 23-4 0 

1 0-221 0-218 0-382 1-68 0-992 | 141 0-48 

2 0-200 0-188 0-247 1-42 0-992 | 121 0-74 

3 0-198 0-196 0235 188 0-992 085 

5 0-190 0-188 0-179 1-05 0-992 8-76 1-33 

6 0-208 0-200 0-179 0-99 0-991 827 146 | 

7 0-208 0-202 0-183 101 0-991 | 8-34 142 | 
10 0-206 0-204 0-190 1-08 0-991 1387 | 
20 0-214 0-212 0-199 1-04 0-990 | 8-55 1-33 
0-216 0214 0-204 1-06 0-968 742 1-82 
80 0-206 0-205 0-211 114 0-978 «152 1-18 

volume : 9-5 mm® + 0-1, equivalent radius : 1:32 mm + 0-01 

0 0-200 0-199 0-822 4:58 0-996 29-0 0 

1 0-186 0-188 0-534 3-22 0-996 19-2 0-29 

2 0-192 0190 | 0423 2-46 0-996 | 14:3 0-59 

5 0-206 0-205 0-314 1-70 0-996 9-55 1-16 

10 0-181 0-178 0-262 1625 0-096 1-24 
25 0-184 0-183 0-270 1-64 0-905 | 9-20 1-28 
50 0-177 0-175 0-271 1-715 0-993 9-05 1-15 
w 0-181 0-177 0-263 1-635 0-998 | 8-63 1-238 
80 0-205 0-314 1-70 0-984 8-20 1-18 


6-95 
6-07 
8-75 
2-54 
2-48 
2-44 
2-54 
2:56 


+ 0-2, equivalent radius : 


0-999 
0-999 
0-999 
0-998 
0-998 
0-997 
0-995 
0-988 


SOLE [7] have found that the partition coefficient 
for iodine between water and carbon tetrachloride 
is independent of the iodine concentration up 
to 0-248 g per 1,000 g of aqueous phase. They 


aq P 
Gwywxne [8] have reported the value of 89-9 at 


in the present calculations. 


expressed 


this temperature, and this value has been used 


The overall mass-transfer coefficients can be 
either on the basis of the distance 
travelled by the droplets, or on the basis of the 
time of contact : 


dN = K pny (C*x — Cg) A - dh == 
= (C*g — Cg): A dt (1) 


Go 
| 
- 108 cm/sec 
radius: 
0 
0-028 
0-040 
0-047 
0-071 
0-078 
0-077 
0-074 
0-074 
0-092 
OL. 
5 
956 
0 
0-018 
0-085 
0-070 
0-072 
0-072 
0-076 
0-081 
0-086 
| Droplet volume : 42 mm* 1-00 mm + 0-02 
| 0188 1-160 1-751 | 20-1 0 0 
0205 0204 1117 | 1522-260 | - - 
16-0 | 0-192 0190 0-648 0027 «149 0-083 
| 160 0-200 | 0-198 0-454 | 0625 100 103 0-066 
157 | 0202 0200 0-468 | | 0610 958 | 107 | 0-070 
| 153 | 0203 | 0202 0-445 | 0-600 918 110 | 0075 
| 148 0-189 | 0186 0-428 | 922 1080-074 
148 O211 | 0-486 | 0-624 892 1038 0.078 
| 
solvent phase mg iodine per ml of 
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CONCENTRATION OF SODIUM OLEYL -p-ANISIOINESULPHONATE ¢/!OOmi AQUEOUS PHASE 


Fic. 2. Overall extraction coefficients on « height and solvent basis, plotted against concentration of surface- 
active agent for 3 droplet sizes. 


or integrated : Key = Kea, (4) 


where v is the droplet velocity. The droplets are 
assumed to have the form of an oblate ellipsoid, 
and it is convenient to express the ratio of 
droplet volume to surface by means of the radius 
of a sphere of volume V, multiplied by a correction 
factor which allows for the departure from the 
spherical form : 


+— 


4 > > 
CONCENTRATION OF SODIUM OLEYL-p~ ANISIDINESULPHONATE, ¢/1OOm! AQUEOUS PHASE 


Fie. 3. Overall extraction coefficients on a time and solvent basis, plotted against concentration of surface- 
active agent for 8 droplet sizes. 
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$10; 
23 mm* | 


10.10) 
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CONCENTRATION OF SODIUM OLEYL~p~ ANISIDINESULPHONATE, 9/100 ml AQUEOUS PHASE 


Fic. 4. Interfacial resistances on a height and solvent basis, plotted against concentration of surface-active agent 
for 8 droplet sizes. 


(5) factors given in column 7 in Table 1 differ from 
unity by less than 1% in most cases. Correction 
for droplet deformation has been reported also 
by Heertyes, Hotve and Tacsma [11)}. 

In the presence of the surface-active agent the 
rate of mass transfer is reduced considerably. It 
(7) is mathematically convenient to define an inter- 

, b\? facial resistance which is additional to the 


(6) 
956 


diffusional resistances in the two liquid phases. 
(:) On an extract basis the equations are : 
a 


dN = = (mC p; -A-dh 
where a and 6 are the long and short axis of the «a) 
ellipsoid. The droplets were photographed and a Cc 
and 6 measured [9]. The resulting correction "ere ada 


| 


CONCENTRATION OF SODIUM OLEYL- p- ANISIDINESULPHONATE AQUEOUS PHASE 


6 710" 


Fico. 5. Interfacial resistances on a time and solvent basis, plotted against concentration of surface-active agent for 
8 droplet sizes. 
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be 
Kea 
OL. 
where 
j= 
(8) 
110", 

| 

OROPLET SIZE 

Ag 

= 

©0223: 

5 


K. P. Linptanp and 8S. G. Tersesen 


The interfacial resistance as defined by equation 
(8) includes any interfacial resistance that may be 
present in the absence of the surface-active agent, 
whereas the resistance calculated from (9), below, 
as a difference between the overall resistances. 
measured with and without the addition, 
expresses the contribution made by the surface- 
active agent only. 


K 


R 


kK without) 
The overall coetticients as well as the interfacial 
resistances ure given in Table 1 both on « height 
and time basis, and are plotted against con- 
centration of surface-active agent in Figs. 2. 3. 
4 and 5. In Fig. 6 the velocity of the falling 
droplets is plotted, also against concentration of 
surface-active agent. 

The reason for using the distance travelled by 
the droplets as a basis for the results. as well as 
the actual extraction time, is the somewhat 
irregular behaviour of the latter observed at high 
concentrations of surface-active agent. As 
pointed out by Cottson and Skixner [10) 
however, the actual time of fall must be expected 
to be more important than the height. 

It will be noted that the coefficients are all 
given on an extract (solvent) basis. and that the 
extraction and resistance coefficients on a raftinate 
basis can be obtained by respectively multiplying 
and dividing by the partition coefficient. Of the 
two liquid films, the aqueous will offer by far the 
greatest resistance to the transfer of sodine due 
to the magnitude of the partition coefficient. 


Disctsstonr 

The results show that an addition of surface-active 
agent of 10° g 100 ml of the aqueous phase has 
reduced the overall extraction coefficient to a 
value which is only 32 to 88%, of the original. 
The percentage reduction appears to be practic- 
ally independent of droplet size within the range 
investigated. Further additions of surface-active 
agent exerted only a very small effect on the 
extraction coefficient, giving a slight increase if 
calculated on a height basis, and a slight decrease 
on a time basis. Up to about 6.10% ¢ 100 ml 
there is a rapid fall. 


The results obtained are similar to those of 
Garner and Hace [8], who found that 1.5 x 10° 
ml 100 ml of the surface-active agent Teepol 
reduced the rate of extraction of diethylamine 
from toluene drops by water to 45% of the 
original value. The present results differ from 
those of GARNER and Hate [3] in three respects: 
the maximum reduction in the rate of extraction 
occurred with a concentration of surface-active 
agent which is about 250 times smaller than theirs, 
the reduction was 67 to 68% against their 55%, 
and further additions of surface-active agent 
were practically without effect, whereas they 
found a marked increase in the rate of extraction. 

In the present investigation it has been possible 
to analyse the results in terms of an interfacial 
resistance coefficient R,. which is given in Fig. 4 
on a height basis and in Fig. 5 on a time basis. 
For not too high concentrations of surface-active 
agent, a single curve was obtained when the 
resistance coefficient based on the extraction time 
was plotted in Fig. 5 against the concentration of 
the surface active agent for the three droplet 
sizes 4-2, 9-5 and 22-8 This result indicates 
that the interfacial resistance is independent of 
droplet size within the range investigated. It can 
be seen from Fig. 5 that at about 6 « 10° g ‘100 
ml the interfacial resistance reaches a constant 
value of 7,400 sec cm (82 sec cm), the figure in 
brackets indicating the corresponding value on « 
raffinate, i.c., aqueous basis. This compares with 
an overall resistance, without addition of surface- 
active agent, of 10,600 (118) and 11,900 (182) 
sec cm for the smallest and largest drops respect- 
ively. The interfacial resistance thus amounts 
to 69 and 62% of the resistance without the 
addition for the smallest and largest droplets. 


A theory of interfacial resistance 

Garner and Hare [8] have shown that the 
maximum reduction in the rate of extraction 
occurred at a concentration of surface-active 
agent which produced only a comparatively small 
reduction in the interfacial tension. They con- 
cluded that there was no direct connection 
between the two phenomena. For the present 
system, Me.uus and Tersesen [9] have also 
failed to find any direct correspondence between 
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interfacial resistance and reduction in surface 
tension. 

The increase in the resistance to mass transfer 
on the addition of the surface-active agent 
occurred at very low concentrations, and it is 
inconceivable that this effect can be caused by 
anything but a layer of surface-active molecules 
adsorbed at the interface. The sharp rise in 
interfacial resistance followed by an almost per- 
fect levelling off, as shown in Fig. 5, strongly 
suggests that the interfacial layer of surface-active 
molecules was almost completed when the bulk 
concentration reached about 6 x 10 g/100 ml. 
This conclusion is at variance with that of 
GaRNER and Hare [8], who state that the 
resistance to transfer reached a maximum when 
the interfacial concentration was far less than 
required for saturation. 

A theory of interfacial resistance will be 
developed on the basis of the assumption that the 
resistance is proportional to the interfacial 
concentration of surface-active molecules. This 
point of view differs radically from that adopted 
by Hutcrinson [2], and by Garner and Hate 
(8). 

The following assumptions are made with 
regard to the adsorption of surface-active mole- 
cules at the interface : 


1. The adsorbed film is of the “ vapour - ex- 
panded " type with strong cohesive forces. 


The rate of adsorption on the free interfacial 
area is proportional to that area and to the 
bulk concentration : 


k, C,(1 — Ap) 


Some of the surface-active molecules striking 
the covered areas are also adsorbed, and they 
find a place at the interface through adjust- 
ments in the film. The accommodation co- 
efficient for the covered areas depends on the 
ease of this adjustment, and decreases rapidly 
with decreasing free interfacial area. This 
leads to the following expression for this 
contribution to the rate of adsorption : 


— Ao)" 
Desorption can only take place at the covered 


areas, and the rate increases rapidly with the 
tension of the film. For an expanded film 
the tension increases slowly with surface 
cocentration when the free area is large, but 
very rapidly when the free area approaches 
zero. This leads to the following expression 
for the rate of desorption : 


k, (1 — 


The concept of * vapour-expanded "’ films has 
been dealt with by Apam [12], whereas MoILuer 
and Co..re [18] have stated that surface-active 
agents, such as efficient emulsifying agents. 
usually form films of a somewhat condensed 
nature. 

The expressions given above for the rate of 
adsorption and desorption lead to the following 
adsorption isotherm : 


_ ke Aol” 
(1— A Af 


Cc (10) 
It will be seen that this equation gives the Lang- 
muir isotherm for n = «2 and m = 0. To fit the 
data appertaining to expanded films, the values 
of n and m must be small, and the equation can 
be simplified to 
k, Af 

It has not been proved that a true adsorption 
equilibrium is always established on the surface 
of a falling droplet, and it is possible that the 
isotherm describes a stationary state rather than 
an equilibrium. 

Assuming that the interfacial resistance is 
proportional to the covered fraction of the 
interface : 


C, = (11) 


Ri = (Aol) Rin (12) 


the following relation between interfacial resist- 
ance and bulk concentration of surface-active 
agent results : 


(13) 


1 — 


In this equation, R',,, is the limiting value of 
the interfacial resistance, and equal to 7-4 « 10° 
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sec,cm, and K the initial slope of the curve in 
Fig. 5, and equal to 5.4 x 10°. The exponent 
was determined from a logarithmic plot of 


against K and found to be 0-1, 
ut) 

The final equation : 


C, = 54 x lo” 


is shown as the fully drawn curve in Fig. 5, and 
is seen to fit the experimental data very well. 


Interfacial resistance and the transfer of turbulence 
Lewis [14] has recently reported an extensive 
study of individual coefficients with a new type 
of transfer cell with agitation in both liquid 
phases. He found that transfer took place by 
eddy, not molecular diffusion, and that the trans- 
fer coefficients for all phases except furfural could 
ky 


be correlated on one curve by plotting — 
v 


1 
against + The turbulence on 


each side of the interface contributes to the 
coefficient, whereas the Schmidt number, which 
is characteristic of molecular diffusion, does not 
enter at all. Lewis found that, while mobile 
films adsorbed at the interface had little effect, 
rigid protein films caused a retardation of trans- 
fer. He explained this as a result of suppression 
of the transfer of turbulence from one phase to 
the other, so that the term in the correlation 
involving the Reynolds number of the other phase 
drops out. 

This explanation of the effect of interfacial 
films does not appear to apply to the present 
results. GARNER and SKELLAND [15] have shown 
that there is no internal circulation in droplets 
of carbon tetrachloride which fall through water 
when the Reynolds number is of the order 
encountered in the present investigation. But if 
there is no form of motion within the droplet, 
there can be no basis for the idea that the ab- 
sorbed film suppresses the transfer of turbulence 
from one phase to another. Further evidence is 
obtained by calculating the droplet side-resistance 
on the assumption that the interior is stagnant. 


The curves given by Heertses, Hoive and 
Tatsma [11], based on the equation of NEwMAN 
[16], were used, and the diffusion coefficient for 
iodine in carbon tetrachloride was estimated to 
be 1-4 x 10% em?/sec by the method of WiLke 
[17]. The resistance was calculated to be 0-6 to 
0-7 sec/em as compared with an observed inter- 
facial resistance of 8 to 9 x 10°sec/cm. This 
result shows that even if complete transfer of 
turbulence to the droplets took place in the 
absence of surface-active agent, and the addition 
inhibited completely this transfer, the resulting 
reduction in the rate of mass-transfer would only 
be a small fraction of that actually observed. 

These considerations show that the observed 
effect cannot be explained by assuming that the 
solvent side mass-transfer coefficient is reduced 
due to suppression of transfer of turbulence from 
the aqueous to the solvent phase. They do not, 
however, rule out the possibility that the adsorbed 
surface-active agent can give to the interface a 
certain rigidity, which could suppress the develop- 
ment of turbulence in the aqueous phase close to 
the interface, and thus reduce the mass-transfer 
coefficient on the aqueous side. On the other 
hand, the observation of Garner and Hate [8) 
that surface-active agents can reduce the rate of 
mass-transfer, even if there is no original internal 
circulation in the droplets, indicates that there 
is no appreciable transfer of momentum across 
the interface, even in the absence of the surface- 
active agent. The evidence on this point is still 
inconclusive. 

The spontaneous interfacial turbulence des- 
cribed by Lewis and Pratt [18] still remains to 
be considered. When the transfer of material 
from one liquid phase to another was accompanied 
by evolution of heat, they observed vigorous 
turbulence and pulsations at the interface. This 
interfacial turbulence would undoubtedly greatly 
increase the rate of mass transfer, and it is 
conceivable that an adsorbed interfacial film 
might reduce or eliminate this turbulence. Accord- 
ing to Davies and Gwynne [8] however, the 
transfer of iodine from water to carbon tetra- 
chloride is endothermic, and in such cases inter- 
facial turbulence has not been reported. This 
explanation therefore does not seem to apply. 
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The nature of the interfacial resistance 


The preceding discussion indicates that the 
interfacial resistance may not be caused primarily 
by changes in the hydrodynamical conditions 
such as transfer of turbulence or interfacial 
turbulence. Some kind of interaction seems to 
occur between the adsorbed film of surface-active 
agent and the diffusing molecules. This inter- 
action could be purely mechanical, or more 
specific, involving forces of a chemical or physical 
nature. Hutcurnson [2] and Garner and Hae 
[8] favour the latter hypothesis, but further work 
is required before a final conclusion can be 
drawn. A comparative study of the action of the 
same surface-active agent on different diffusing 
solutes is in progress in this laboratory. The 
present system is particularly suitable for such 
studies, because the interfacial resistance reaches 
a practically constant value after the addition 
of very small quantities of the surface-active 
agent. Although the mechanism of the retardation 
phenomenon is not yet known, the results of the 
present work show that the effect can be described 
in terms of an interfacial resistance proportional 
to the interfacia® concentration of the surface- 
active agent. 


Velocity of fall for 3 droplet sizes. plotted against concentration of surface-active agent. 


Retardation of the rate of droplet fall 


The retardation of the falling droplets caused by 
the addition of the surface-active agent is shown in 
Fig. 6, and is analogous to that reported by 
Stuke [6] for rising gas bubbles. The sharp 
primary retardation appears to reach its maximum 
value at a concentration of surface-active agent 
about one-third to one-half of that for the inter- 
facial resistance. The present results throw no 
light on the nature of this effect, and no explana- 
tions can be offered in addition to that already 
put forward by Stuke. 
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NOTATION 


A = Interfacial urea cm? 
Ag = Effective cross-sectional area of adsorbed 
surface active molecules em? /g mol 


b/a = Ratio of short to long axis of ellipsoid 
Cg = Concentration of iodine in solvent phase mg /cm* 


Cg* = Equilibrium concentration in solvent phase 
3 


mg/em 
Cp = Concentration of iodine in aqueous phase 
mg/em® 
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Cr’ = Original concentration in aqueous phase mg/em® N ~ Amount extracted 
Average of initial and final Np, = Reynolds number 
aqueous phase mg/cm 100 
‘s ~ Concentration of surface-active agent in aqueous p = Approach to equilibrium — ‘Cre 
phase 8/100 ml Ry, = Interfacial resistance on a solvent and height 
J = Correction factor given by equation (7) basis dimensionless 
h ~ Average effective column height om Ry», = Interfacial resistance on a solvent and time basix 
Kya) ~ Overall transfer coefficient on a solvent and sec /em 


height basis dimensionless 1 
= Limiting value of interfacial resistance sec 
Kp Overall transfer coefficient on a solvent and time r = Equivalent radius of droplet em 


basis cm/sec 
k, = Constant V = Droplet volume em 


k, = Constant v = Droplet velocity em /sec 
m — Distribution coefficient : I’ = Surface concentration of surface active agent 

mg iodine per mi solvent phase g mol/em? 

mg iodine per ml aqueous phase » ~ Absolute viscosity poise 

» = Kinematic viscosity cm? ‘sec 
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Forced heat convection in cylindrical channels: 
Some problems involving potential and parabolic velocity distribution 
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Abstract—Analytical solutions are presented for two kinds of problem in the convection of 
heat by a fluid. The first problem is that of the uniform and constant generation of heat within 


the fluid; the second is that of temperature equalization when the inlet temperature has one 
constant value for part of the radius and another value for the rest of the radius. Steady-state 
solutions are developed for both plug flow and laminar flow ; the wall of the tube is assumed 


to be isothermal. 


Résumé—L ‘auteur présente des solutions analytiques de deux sortes de problémes sur la con- 
vection thermique par un fluide: le premier relatif & la production constante et uniforme de 


chaleur au sein du fluide ; le second concernant |'égalisation des températures quand la tem- 
pérature interne a une valeur constante pour une partie du rayon et une autre valeur pour le 
reste. Il donne des solutions 4 l'état stationnaire soit pour l’écoulement en piston soit pour 


THE merit of the analytical approach to problems 
in convective heat transfer has been demon- 
strated in the sixty-five years that have elapsed 
since the pioneer theoretical studies by Gratz 
and by L. Lorenz. The design of equipment, 
the prediction of its performance, and the planning 
and interpretation of experimental research are 
often conducted more effectively when the 
implications of the heat equation in the particular 
problem are clearly stated. 

Two different kinds of problem in forced con- 
vection are solved here. In both of them, a fluid 
flows through a cylindrical tube. The tube wall 
is isothermal, the properties of the fluid are 
independent of temperature, and we assume the 
fluid velocity to be entirely in the axial direction 
and to be steady. In the first problem, the fluid 
inlet temperature is uniform and there is a con- 
stant heat generation per unit volume; in the 
second, the inlet temperature is uniform at one 
value for part of the radius and at another value 
for the rest of the radius. Both constant and 
parabolic velocity distributions are considered in 
each of the problems. The steady-state local 


I'écoulement laminaire, en supposant que la paroi du tube est isotherme. 


* Present Address ; Esso Research and Engineering Co., P.O. Box 121, Linden, New Jersey. 
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temperatures are calculated and presented 
analytically and also graphically. The solution 
of the first problem should be helpful in estimating 
temperatures developed in chemical and nuclear 
reactors, and that of the second problem in 
estimating the average life of non-uniformities 
of the temperature in such streams, 


I. Unirorm Heat GENERATION 


The differential equation of heat conduction in a 
moving medium is 


=q/pe+aA*T. (1) 


Assuming symmetry about the axis of the 
cylinder, neglecting axial conduction of heat 


FF 
(assuming <= at the steady- 
state (1) becomes 
oT YT 
(2) 


The fluid enters at z = 0, at the uniform tem- 


5 
| 
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perature Ty, and the wall (r = 8) is at the constant 

temperature T,. Define @as T T,. Then the 

boundary conditions are : 
atr=0,0=0,=T,—T, (3) 
atr=s, @=0. (4) 


A. Potential Flow 
Here }, is a constant, V, and (2) becomes 


M0 +a 


(5) 


The homogeneous equation (4 = 0) has the 
solution [references 1, 2] 
ae g2 @)\ (6) 
B,J (B,) Jo (But) exp ( sv 
The values of 8, are the positive roots of 
J(8,) = 0, w =(r/s); J, and J, are Bessel 
functions of the first kind and zero and first 
order respectively. 

The solution of (5) has the form 


Nn JJ exp ( a? *) 


This automatically satisfies (4); N, must be 
computed from boundary condition (8) : 


Multiply (8) by J,(8,,) wdw where need 
not be the same eigen value as 8,, and integrate 
from w= 0 to w= 1. 

For each n, we have the result 


1 1 
0 { dee — N, Jy(Bqte) de 


1 1 
As* 

The following properties of J, are introduced in 

(9) to produce (14) : 


1 


1 


1 
B, 


As* J 


and 
and finally, from (7), 
x 2 A? 1 


exp (— + (16) 


The table below is for use in computing the first 
three terms of (16). 


n Bn J(B,) 
1 240 0-520 
2 5.52 — 0-840 
8 8-65 0.272 


Fig. 1 is a plot of axial temperature (w = 0), 
from equation (16), using appropriate dimension- 
less groups. The average (or bulk) temperature 
at any axial station is another important quan- 
tity. It is defined as 


fears V, dw 
6. = (17) 
V, dw 


Fig. 2 expresses the difference between the 
axial temperature and the mean temperature. 


B. Laminar Flow 
In laminar flow, V, is a function of w : 
V, = 2¥,,(1 —w*) 
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Axial temperature in a tube. Potential flow, 
uniform heat generation, isothermal wall. 


Fic. 1. 


As? As? As* 
1) =0; (2) =1; (8) =4; — = 6. 


The differential equation here is 


, #1290 
2 5 ("+2 w dw 
with the boundary conditions of eqs. (3) and (4). 
The homogeneous form of eq. (19) was solved by 
Graetz [2] discussed in the modern literature, 
{1}, [8], [4]), who also presented equation (6). 

Graetz found the solution of the homogeneous 


+4 (19) 


— Ba) 


R (w, B,) = Bzn (w 


(21) 


where B, = l, B, = Ban 
and the 8, are the roots of R(1, B,) 


Fic. 2. Difference between axial and mean temperature. 
Potential flow, uniform heat generation, isothermal wall: 


ay) =0; (2) = 4; (s) = 6. 
aby 


=0. This function was studied by Graerz, 
and a table of its numerical values follows. 
The first three zeros are (4): 8, = 2-704, 8, = 
6-68, B, = 10-67. 

A solution of eq. (19) is 


a= NR(w, exp (— 


+ As —w*) (22) 
ta 


The N, of eq. (22) are calculated from eq. (3), 
since boundary condition (4) is automatically 
satisfied by (22). At z= 0, 


0 =0,= EN, Rw, + w*) (23) 


n 
Multiply eq. (28) by R,, w (1 — w*) dw, where m 
need not be equal to n. By R,, or R,, we will 
mean R (w, §,,) or R(w, B,) respectively. Inte- 
grate from w = 0 tow = 1, 
For each n, 


'5 
| | 
| 
6 
O4 
| 
0 0 
O3 4 
+—— 
a! 
+ - + 
OL. ke re) O2 “O06 O8 
5 
956 
— 
@ 
— 
vB 
«@ *) 
con 20 
exp( 2sVis (20) 
= 
n=O 
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(1 — dw 


= R,, w w*) dw 
(24) 


« 
w*) 


| (1 
These properties of R are used in (24) : 


dw 


R,w(l = 0, men 


(25) 


R260 w*) du ae 


i 
Oo 


= 
SV 


| — w*) de 


| R,w* (1 — w*) dw 


2 { 
+ R, (28) 
Bn 


Equations (25) to (28) are due to GRarTz ; 
follows from an integration by parts. 
These lead to 


(28) 


w* dR, 


The first terms of the series solution of (22) are 
then 


(29) 


T 
LN 


o8 10 


+ 
| 
| 


oo O2 


Fic. 3. Axial temperature in a tube. Laminar flow, uniform Fic. 4. Difference between axial and mean temperature. 


heat generation, isothermal! wall: 
aby 


A? 
(4)-—— = 6. 


Laminar flow, uniform heat generation, isothermal wall: 


A? A* As? 
(1) —— (2)— =4; (8) — = 6. 
aby 


1 (dR, 
v | 
1 
1 
wel 
1 
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— 2202 40" ) ~3.66 *) 
sV 8 


+ (0.588 7; (30) 
Fig. 8, is a plot of axial temperature (w = 0), 
from equation (30), using appropriate dimension- 
less groups. Fig. 4 expresses the difference 
between the axial temperature and the mean 
temperature (eq. (17)). 


2. Temperature Mixine Between 
CONCENTRIC STREAMS 

The rate of * mixing-out " of non-uniformities of 
temperature or of composition in a flowing stream 
is often important to the analysis of processes. 
The calculations that follow are equivalent to the 
rate of equalization of temperature due to the 
diffusion of heat alone. In most practical] situa- 
tions, the initial non-uniformity of temperature 
induces non-uniform velocities that expedite the 
temperature equalization. For this reason, the 
actual rate of temperature mixing should always 
be at least as great as indicated here. 

Consider a stream of fluid entering a cylindrica! 
tube of radius S. At z = 0, the stream tempera- 
ture is T, from the axis (r = 0) tor = a(s >a); 
itis T,fromr = ator = s. The wall temperature 
remains constant at 7',, and the assumptions of 
eq. (1) apply here as well. 

The differential equation is 


_ 190 (at) 
S? 


with the boundary conditions 


Atx = 0,0 = 0, = T, — T,. (a;s) >« >0 
T, —T,, 1< w < (a/s) (82) 


Atw=10=0 (33) 


A. Potential Flow 
JV, is a constant. and (31) becomes 


Forced heat convection in cylindrical channels 


which has the solution 


sls 

where the 8, are the positive roots of J, (8,,) = 0. 

The N, are determined from eq. (32). Subject 

to certain restrictions on f(w) that are satisfied 


here, when 


N, J, (8, w) (36) 


2 
N, [J (B f (w) J, (B,, dw (37) 
1 n 


f(w) = 


Here 
/s 


» a 
N. = — - yd 


0, | w) (38) 


a/s 


2 a a 
(39) 
And finally, 
0 2 —% a 
J 
6, B, (B,)}? F (“ A, (2, ‘) 


(2 J, (8, exp — 4 (40) 
Fig. 5 indicates the approach of the temperature 
of the fluid to that of the isothermal wall when 
(a s) is 0-5, and (°: = oy is 0-5. The radial 


temperature distribution is plotted at x = 0 and 
at two values of the dimensionless distance down- 


a 
stream, ‘). 
(. Ves 


B. Laminar Flow 
The problem is 


0 12 
dz 3 (23+. (34) 
5 
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1 
|| 
17 


Table 1. 


Graetz's R, Function (4) 


Ww Ro R, Ry 
00 1 1 1 

0-982 0-892 0-726 
0-2 0-929 0-605 0-153 
0-38 0-846 0-234 ~ 0815 
0-788 0110 — 0-892 
0-5 0-615 — 0-840 — 0-142 
0-6 0-483 — 0-482 0-170 
o-7 0-851 — 0-398 0-382 
0-8 0-224 — 0-284 0-308 
0-9 0-107 0-141 0-168 
10 0 0 0 


with boundary conditions of eqs. (82) and (88). 
The solution is 


2 a @ 


n (42) 


N, R (w, B,) exp ( 


where R(1,8,) = 0. The N, are from eq. (82) : 


= N,R,(w) (48) 


Multiply eq. (48) by R, w(l — w*) dw, where 
m need not be equal to n. Integrate from w = 0 
tow = 1. For each n, using eq. (25) to (28), 


ass 


+6 —w*)dw (44) 


2 


a (dR, /dw)..4 


+ (45) 


() is most conveniently calculated by 
dw 
graphical differentiation of R,, or else by a 


graphical integration : 


a/s 


de wea/s a/ 
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The first three terms of the solution to eq. (85) 
are 


dR, 


dw 
1-01 


Ry (w) exp (- 3.66 


s\ O 1-33 (47) 


R, (w) exp (— 228 

a 0, —9 (dR,/dw), 

az 


Fig. 6 expresses the temperature equalisation in 
laminar flow in a tube having an isothermal when 
(a/s) is 0-5 and (*) is 0-5 (the same con- 


ditions as for Fig. 5). 


az az ae 
—- = 00; —-=01; (8)—- = 05. 
Os 
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5 
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| 
lw | 
— 
| | 
0-4 — | 
02 
= 
| 

Ww 
ee Fio. 5. Equalization of temperature. Potential flow, 

(a/s) = 05, | ———]} = 
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NOTATION 


A = 4q/pec 
= Bessel functions of the first kind and 
zero and first order 
coefficient of characteristic function in 
various equations 
= Graetz’s function, defined by eq. (21) 
= temperature, initial temperature at 
r < a, initial temperature ata < r<s, 
initial temperature, wall temperature 
mean velocity, local axial velocity 
radial coordinate at discontinuity of 
initial temperature 
¢ = heat capacity 
q = heat generation per unit volume and 
time 
r = radial coordinate 
8 = tube radius 


w = (r/s) 
9 Kaualizati tin wv = axial coordinate 
56 “qualization perature. Laminar flow, a = thermal diffusivity 


(a/s) = 0-5, = 05. 8, = a zero of Jy or of R,, 
% p= density 
7 = time 
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The simultaneous transfer of heat and mass between water and moist 
coal gas 


G. S. Crrss and T. Ne.son 
Fulham Research Laboratory, North Thames Gas Board 


Abetract—The Graphical Method of Micx.ey (8) for the counter current contacting of moist 
air with water is modified to apply to the system moist coal gas- water, over the range of saturation 
temperatures 82°F to 180°F. Over this range the Lewis Number Ag/kgs cannot be considered 
constant and the log mean partial pressure of the non-diffusing gas varies widely from unity. 
A driving force is derived, which incorporates these variables for use in conjunction with the 
ges film mass transfer coefficient. This necessitates the use of two equilibrium lines as in the 
method of Lewis and Wars [5). 

Although no suitable data are available for truly counter current flow for the system considered, 
the method has been applied to test data on a cross flow water-tube condenser. The results 
are compared with those obtained by the well known method of Cotsurn and Hovoen [1). 


Réesumé— Les auteurs ont transformé la méthode graphique de Mickey valable pour l'échange 
entre de l'eau et de l'air humide circulant & contre-courant afin de l'appliquer au systéme 
eau-taz de houille, au-dessus de lintervalle des températures de saturation (82°F-180°F). 
Au-dela de ce domaine, on ne peut pas considérer le nombre de Lewis hg/KgS comme constant et 
la moyenne logarithmique de la pression partielle du gaz diffusant s'écarte largement de l'unité. 

Les auteurs déduisent une équation de la force motrice qui comprend variables et en méme temps 
tient compte du coefficient de transfert massique du film gazeux. Cela nécessite l'emploi de 
deux lignes d'équilibre comme dans la méthode de Lewis et Wurre. 

Quoique dans le systéme envisagé, il ne soit pas possible d'avoir des données valables pour 
un écoulement réellement a contre-courant, ils ont vérifié par cette méthode les valeurs relatives 
& un condenseur & tube d'eau avec circulation externe & angle droit. Ils ont comparé leurs 
résultats aux résultats obtenus par la méthode bien connue de Cotpunn et Houcen. 


The method may also be employed to deter- 
mine values for individual film coefficients from 
performance data on actual plant. 


1. INTRODUCTION 


Since Cotsurn and Hovcen [1] first published 
their paper on the ‘ Design of Cooler Con- 
densers for Mixtures of Vapours and Non- 
Condensing Gases,’’ several less rigorous but 2. Tue DEVELOPMENT OF A 
simpler methods have been proposed. Co-BuRN GraPnicaL METHOD FOR THE 
[2] in a recent paper stressed the need for such DesioN oF Piant INVOLVING 
methods, based preferably on the terminal THE SIMULTANEOUS TRANSFER 
conditions only, and suggested that the most or Heat anv Mass 


hopeful direction for investigation was the 
method proposed by Micx ey [8]. 

By extending the method of Micx.ey to a 
system other than air-water having high vapour 
concentrations, a graphical method is developed 
which is as rigorous as that of Co.surn and 
Hovucren. Moreover, the method does not 
assume that the relative rates of heat and mass 
transfer are such that the gas remains saturated 
throughout the apparatus, an assumption which 
is shown to be invalid. 


If in countercurrent operation of contacting 
moist gas with water, a small section of transfer 
area dA is considered; five basic differential 
equations for heat and mass transfer may be 
written down. 

The mass balance relating to the transfer of 
the water between the gas and liquid phase is 
given by 

G dH, = dL (1) 


The complete heat balance over the section, 


__ 
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The simultaneous transfer of heat and mass between water and moist coal gas 


expressed as changes in enthalpy based on gas 
and liquid water at t, (82°F), is given by 


G dig =I dt, + dL Cy (ty to) (2a) 


The enthalpy change of the gas comprises four 
terms 


G dig = Gsdtg + GdHgA + 
G dHg ty (tg — + dHg — ty) (2b) 


where Gsdtg is the sensible heat change in 
the bulk of the gas. 

G dHg A, is the latent heat of the mass of 
water transferred across the interface at 
temperature 

G dHg cy (tg — t,) is the sensible heat change 
in the water vapour in passing between the 
bulk of the gas and the interface. 

G dHg cy, (t; — to) is the sensible heat or 
enthalpy of the water transferred above 32°F. 

The three remaining basic equations are as 
follows :— 

The rate equation for sensible heat transfer 
across the gas film 


— Gadtg = hg (tg — 4) dA (3) 
The rate equation for mass transfer across the 
gas film 

=kg — 4A (4) 


The rate equation for total heat transfer across 
the liquid film 


Ley dt; = hy (ty — 4) dA (5) 


The difference between dlc, (t; in 
eq. (2a) and GdHg c, (t; — to) in eq. (2b) is 
G dHg ¢, (t; —t,) the heat given up by the 
condensate in cooling down to the bulk liquid 
temperature ¢,. This quantity is small com- 
pared with GdH, A, and may be neglected. 
The significance of the quantity G dHg cy 
(tg — t;) is discussed below. 

The total sensible heat flux at a point in the 
gas film is given by 

G dH, 
dA 


where z is the fractional distance through the 
gas film corresponding to temperature. CoL- 


cy (t — &) (6) 


BURN and Drew [4] integrated this expression 
across the gas film to give 


= he (ta — 4) (7) 


where 
kg Cy 
G dHg cy/hgdA 


he 


For moist coal gas upto a saturation temperature 
of 180°F, a lies between 0 and 1, and under 
these conditions —e¢“ approximates to 
1 +0-54a. Thus the contribution to the 
sensible heat flux by the diffusing vapour 
approximates to 054 GdHgcy (tg —t) or 
G dHg (tg — t;)/4. This term is usually very 
small compared with G dH, A; and may safely 
be ignored. 

The heat balance for the total heat transfer 
across the interface is therefore given by 
equation (9), 


G dig = Gs dt + GdHgi,=Le, dt, (9) 


(H;, H,) (8) 


L c, dt, is @ closer approximation to the total 
heat transferred across the interface than G di,. 

In order to avoid laborious trial and error 
calculations employing the above equations, 
methods have been developed for combining 
the heat and mass transfer driving forces across 
the gas film into a single potential difference. 
For the system air-water, a fortuitious relation- 
ship exists over the range of cooling tower 
operation. The ratio hg ‘kgs (the Lewis Number 
b) approximates to unity under these conditions, 
and eqs. (3) and (4) may be combined as follows 


— G(sdtg + dHg) = kg [(stg + He) — 
(st; + Ay dA (10) 


Assuming s is constant, the left-hand side 
of the equation is equal to — G dig, and the 
right-hand side approximates to the difference 
in enthalpy between the moist air and that in 
equilibrium with water at the interface tem- 
perature (t;). This equation forms the basis 
of the graphical method of Mickiey [3] for 
the air-water system. 

For cases where the Lewis Number is not 
equal to unity but may be considered constant, 


re 
= 
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Lewis and Wurre [5] have proposed a * modi- 
fied enthalpy method.” In the gas industry. 
gas coolers handle highly saturated gas over a 
wide range of temperature. Under such con- 
ditions the humid heat, on the basis of unit 
volume of dry gas and consequently the 
corresponding Lewis Number, vary by as much 
as 100%. Any method therefore that assumes 
constant humid heat is not applicable. Further- 
more, humidity as a potential for mass transfer 
is a valid approximation, only when p_ the 
partial pressure of the diffusing vapour is 
small compared with P — p, the partial pressure 
of the non-diffusing gas. 

At 180°F the vapour pressure of water is 
15-29 in. Hg. Thus for saturated gas at 180°F 
and 30in. Hg. total pressure, the partial 
pressure of the water vapour exceeds that of 
the gas, and humidity difference cannot be used 
as an approximate driving force. The expression 
for humidity driving force is H,; — Hg 
E 

P—-p, P—De 
or potential difference for mass transfer, under 
conditions of steady state diffusion of one gas 
through a second stagnant gas, is given by 
Snerwoop and Picroxp [6] as In 

- 


| The true driving force 


This may be written as In ( — in| 


P—p; P—pg 
Comparing the two expressions for driving 
force, the true potential for mass transfer is the 
humidity multiplied by a factor <. 


(11) 
P—p 


In (1 +2) 
and if the right-hand side of the equation is 
expanded and terms of 2* and greater powers 
ignored, then « = 1. Values of « are tabulated 
below for saturation temperatures from 82°F 
to 180°F. 


Equation (11) is of the form x 


32 #0 8 100 120 14 160) 180 
a 0-997 0-996 0-991 0-982 0-969 0-925 0-895 0-819 0-686 


Humidity as a driving force therefore is a 
very good approximation up to saturation 
temperatures of 90°F, but the error at higher 
temperatures may be illustrated as follows. 
Humidity may conveniently be expressed in 
the gas industry as pounds of water per 1,000 
eu. ft. of gas at S.T.P. (60°F and 80 inches 
saturated). To simplify the units of humidity 
and humid heat, this quantity of gas will be 
termed a Standard Gas Volume, abbreviated 
to S.G.V. 

Consider saturated gas at 180 °F in contact 
with water at 160°F, 

From Table 111 in Sprers 5th edn. [7]. humidity 
driving force between saturated gas at 180°F 
and gas in equilibrium with water at 160°F is 
equal to 

49-19 — 22.39 = 26.80 Ibs. 'S.G.V. 


True driving force = (0-686 x 49-19) — (0-819 
x 22-39) = 15-31 Ibs. (S.G.V.). The error is 
therefore 75%. 

The rate equation for mass transfer is more 
correctly expresesd as 


— GdHg — Hj dA (12) 


Thus the mass transfer coefficient k, used in 
conjunction with aH as a driving force, is 
equivalent to the kg Pyy used in conjunction 
with partial pressure as a driving force. The 
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latter form is that normally encountered in 
mass transfer problems [6]. 

The true potential aH is a unique function 
of humidity H and may be added to the 
humidity temperature diagram, more familiarly 
known as a psychrometric chart (Fig. 1). 

Equations (8), (9) and (12) may be combined 
to give the following expression 


— Gdig = [he (te —t) + 
kg —%,H,))dA (18) 


In order to combine the driving forces for 
heat and mass transfer, one further relationship 
is required in addition to the five basic differ- 
ential equations. This is the relationship be- 
tween hg and kg, and it is assumed, with 
theoretical justification, that a constant value 
for the ratio hg/kg may be assumed throughout 
the apparatus. and CoLsurn [8] 
have demonstrated that for a variety of flow 
conditions the ratio hg/kge is substantially 
constant. The heat capacity c is expressed on 
the basis of the total number of moles. The 
molar heat capacities of coal gas and water 
vapour are 7-7 B.Th.U/lb. mole and 8.0 
B.Th.U/lb. mole respectively. Thus_ the 
variation in c is small over the whole range of 
humidity. 


E = hg/kg = a constant (14) 


Equation (18) may therefore be re-written as 
follows : 


— G dig = ke { [E (te — to) + %@ He) — 
[E(t —t) +A dA (15) 


Applying the equation to saturated gas and 
introducing Ag 


— Gdig = kg{ [E (ts — to) + As xs Hs] 
+ (A—As) asHs —[E(t, A (16) 


Thus the driving force is the difference between 
two potentials, plus a small deviation 
(A, — As) as Hs. 

It will be shown later, by reference to an 
example based on actual test data, that this 
deviation is small, and it will be ignored in the 
further development of the method. 
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The form of the uncorrected enthalpy transfer 
potential 


= E(ts —t) + Asas Hs (17) 


is such that if E and Ag are constants, points 
of equal potential are connected by a straight 
line of slope — E/Ag on a plot of aH against 
temperature. With reference to Fig. 1, the 
condition of unsaturated gas may be repre- 
sented by a point (tg, ag Hg), to the right of 
the saturation line on the ordinate through 
(tg, Hg). This point may be related to a 
saturation temperature ts, on the saturation 
line ag Hg vs. ts, by means of a line of slope 
— E/As. The enthalpy transfer potential 0, 
of the unsaturated gas (tg, Hg) is therefore 
equal to @, that of saturated gas at a tem- 
perature fy. 

Humidity determination by means of a wet 
and dry bulb hygrometer involves the ratio of 
the heat and mass transfer coefficients. Ift, and 
tw are the dry and wet bulb temperatures, as 
recorded on the hygrometer, then the humidity 
H, of the gas is given by the relationship 


he (tg — tw) = ke Aw (Hy — He) (18) 
or 

My 

tw — te Aw 

For the temperature range over which E is 
determined, « is approximately equal to 1, and 
if the wet bulb relationship is to be extended to 
temperatures above 90°F, then (Hy — Hg) 
should be replaced by (a. Hy — ag Hg). Lines 
of ‘ constant transfer potential’ are therefore 
lines of ‘constant wet bulb temperature’ if 
the value of E is the same for both the apparatus 
under consideration and the wet and dry bulb 
hygrometer. In packed towers where, if wetting 
is incomplete, the area available for heat trans- 
fer exceeds that available for mass transfer, 
then E varies if the coefficients are expressed 
on the basis of unit volume of packing. How- 
ever, at high liquor rates, Hensex and Trey- 
BAL [9] found that the ratio of coefficients 
agreed with the psychrometric ratio for air- 
water used in conjunction with the wet and 
dry bulb thermometer. The rate equation for 


total heat transfer in the gas film may therefore 
be written as 


~ G dig = kg (0g —0,) aA (20) 
or 
Gf dic 


fer 
ke J (0g — 4) 


(21) 
The two remaining relationships required for 
a graphical solution of problems involving the 
simultaneous transfer of heat and mass are 
given by the combination of equation (20) 
with the other rate equations. 

Combining (5), (9) and (20) 


th 


am (22) 
(t, — ti) ke 

Combining (3) and (20) 
dig kgs (0¢ (28) 


he (te — 4) 
The significance of these equations will be 
explained with reference to Fig. 2, a plot of 
the enthalpy of moist gas at saturation (ig) and 
the transfer potential at saturation (@,) against 
the saturation temperature. Line AB is the 
operating line representing equation (2) and 
having a slope of L.G. The line of constant 
enthalpy or isenthalp AC represents the en- 
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thalpy of moist gas at the hot end of the 
apparatus. 

Points A and C correspond to the water and 

gas temperatures respectively. The line of 
constant potential A’C’ represents the transfer 
potential corresponding to the terminal con- 
dition of the moist gas. (f,,), and thus @¢, Is 
obtained from (ig,, tg) by a “ constant wet 
bulb temperature line“ on the pyschrometric 
chart. . 
The point of intersection D, of a line of 
slope — h, kg through A’, with the potential 
saturation line, gives the conditions at the 
interface. 

Line C’D has a slope of (@g — 8,) (tg — t)- 
Equation (20) may be re-arranged as follows 

b dt; (24) 

dig (Ag — 

If the isenthalp FG is the enthalpy of the gas 
after a small change dig then the intercept GJ 
cut off by the intersection of a line parallel to 
C’D with FG, is equivalent to 1/6 dt. Thus 
the correct gas temperature, after the first 
increment of transfer, is given by point E 
where EG = bJG. Line CE is the instant- 
aneous direction of the gas path and the 
required value of } is that corresponding to 
(Yeas 


CONSTANT 
ENTHALPY LINE 


HUMIDITY, 


| ASS 


The representation of enthalpy on the psy- 
chrometric chart by a straight line relationship 
is an approximation, and the percentage 
deviation from the true value may also be 
plotted on the chart as a series of curves. From 
a knowledge therefore of the enthalpy and the 
temperature of the gas; the humidity, and 
thus the transfer potential, can be obtained as 
shown in Fig. 8. This diagram illustrates the 
connection between the saturation temperature 
on a constant enthalpy basis (ts); and the wet 
bulb temperature of the gas (ts),. This con- 
version has to be carried out for each increment 
of enthalpy. For saturated gas the two 
saturation temperatures are of course identical. 

The procedure for the graphical method is 
therefore as follows : 


1. Determine the overall enthalpy changes in 
the gas and the liquid. 


2. From these data, plot the operating line 4B 
on the enthalpy-temperature chart (Fig. 2). 


8. Plot @, vs. ty on the same chart. 


4. Sub-divide the operating line into ten equa] 
increments (ig;, ten) 


5. Determine the value of (ts,), from (ts,)é 
(Fig. 8) and hence the value of @,, (Fig. 2). 


6. Construct a tie-line of slope — A, /kg and 
locate Oi, ti (Fig. 2). 


7. Join point (Ag), tg,) with point (@,,, 4) and 
construct a line parallel to this line through 
ig), tg), to intersect a line parallel to the 
abscissa with ordinate igs. Locate point tg, 
on this line by the method indicated (Fig. 2). 


8. Repeat steps 5, 6 and 7 for the remaining 
increments. 


9. Tabulate the values of and apply 


= 
(0g — 9) 
Simpson's Rule to obtain [dig/(@g — 0,) the 
number of individual gas film transfer units 
required. 


The transfer area required is then given by 
eq. (21) 


G di 
Aw as 
— 4) 


The simultaneous transfer of heat and mass between water and moist coal gas 


The reverse procedure to determining A from 
individual coefficients, is to determine the 
coefficients from tests on actual plant of known 
transfer area. The problem then arises as to 
what value to assume for the tie-line slope. 
MICKLEY claims that this value may be found 
from a single test, by adjusting the slope until 
the humidity of the outlet gas agrees with the 
observed value. The disadvantages of this 
method are as follows : 


1. The locus of the gas condition is dependent 
on the method of selecting the increments. 


2. Small errors in the outlet gas condition 
correspond to large variations in tie-line 
slope. 

3. If the gus becomes saturated at any point 
it is almost certain to leave the unit in a 
saturated condition irrespective of tie-line 
slope. 


An alternative method is to eliminate the 
liquid film resistance for a given set of flow 
conditions by adjusting the water temperature 
so that it is constant throughout the apparatus. 
Under these conditions, the loss of sensible heat 
from the gas is balanced by the gain in latent 
heat, and thus the method is only applicable to 
direct contact apparatus such as packed towers. 
Furthermore, the method requires adiabatic 
conditions and the measurement of small 
driving forces. It is in fact suitable only for 
experimental plant and not for checking the 
design of plant under operating conditions. 

A third method is suggested by expressing 
the overall coefficient for enthalpy transfer in 
terms of kg and h,. 


G dig = Kg (Og dA 
= kg (Og = hy (t; tr) dA 


Now (t; —t,) = ~ (0; — 0.) where m is the 


slope of a chord on the saturatign line (@y vs. 
ty) 


G dig ke ht (25) 


Equation (25) is identical with the expression 
for the overall coefficient in gas absorption, 


— 
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where m is the slope of the equilibrium line. 
For the system under consideration, m increases 
rapidly with temperature, thus there is a 
corresponding decrease in the value of Ky. This 
is a particular case of the general axiom, that 
if the tie-line slope is assumed is too steep, the 
value of kg so determined will decrease with 
temperature. If therefore the unit can be 
divided into two sections of equal area, the 
tie-line slope can be adjusted until the calculated 
values for kg in the two halves are equal. For a 
multipass tubular condenser this is possible if 
the interpass water temperatures are recorded. 
The water velocity must however be such that 
thermosyphoning, as described by 
[10] is avoided. In direct contact apparatus, 
where it is not possible to measure intermediate 
temperatures, the humidity of the gas could be 
determined by aspirating a sample through a 
wet and dry bulb hygrometer. This third 
method of ascertaining the correct tie-line 
slope should therefore be more applicable to 
test data on full scale plant obtained under 
working conditions. 


The graphical method which bas _ been 
developed requires a knowledge of the enthalpy 
and the transfer potential of moist gas and the 
relationship between these two properties. This 
relationship is obtained from the psychrometric 
chart, the construction of which is described 
below. 


8. CONSTRUCTION OF A 
PsycuHROMETRIC CHART FOR 
Moist Coat Gas 


The psychrometric chart is based on a plot of 
humidity vs. temperature, and these data may be 
obtained for coal gas from Spiers 5th edn., 
Table 111. A convenient unit for gas industry 
purposes, is pounds of water vapour per 1,000 
cu. ft. of gas at S.T.P., abbreviated to Ibs. /- 
S.G.V. When the humidity saturation line has 
been constructed, the plot of ag H, may be 
added using the appropriate values in Table 1 
for a. 


As stated previously, lines of “ constant 


transfer potential '’ are straight lines on the 
psychrometric chart of slope — E/Ay. 


E is 


Table 1. 


(tg — fo) Ag as Hs Os, 
B.Th.U./ B.1h.U./ B.Th.U, / 
(S.G.V.) | (S.G.V.)| (S.G.V.) 


1,200 
2,011 
8,275 


146 416 

776 
1,053 
1,830 
1,607 
1,884 
2,161 
2,438 
2,715 
2,902 
2,269 
8,546 
8,821 
4,101 


8,050 
8,982 
5,134 
6,547 
8,317 
10,515 
18,340 
16,990 
21,880 
28,420 
87,510 


1,244 
1,427 
1,610 
1,793 
1,976 
2,159 
2,842 
2,525 
2,710 


E, = 27-7 B.Th.U/(S.G.V.) (°F). 


Ey = 18-8 B.Th.U/(S.G.V.) 
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1 2 3 8 | | 12 12 

cy Hs E, | | + 

ts Hs as As (ts ts — &) | O45 Hy 

Ibs. / B.Th.U./ B.Th.U./B.Th.U. //B.Th.U.//B.Th.U. | B.Th.U. //B.Th.U. / 

$.G.V.)| (S.G.V.) | (S.G.V.) | (S.G.V.) (S.G.V.)_ 8.G.V.)| (S.G.V. 
F |G. ) ( )| (S.G.V.) | ( )) (S.G.V.) ( ) 
40 | 0-390 | 0-996) 160 4205 «140/581 638 | | 20-18 
50 | 0573 0-904) 360 617 | 464, 982 1,106; 936 | 20-26 
60 | 0-829 0-901 560 802 | 10-48 | 1,462 1,667 | 1,888 | 20-37 
70 | 1188/0087) 760 1,978 | 21-8 | 2,085 2,294 | 1,926 | 20-53 
so | 1-669 0-082) 960 | 1,796 | 360 | 2,792 2,508 | 20-76 
90 | 2882 0-976) 1,160 | 2,509) 609 8,730 8,436 | 21-05 
100 | 8-234 | 0-060| 1,860 | 8480, 088 | 4,039 4,404 | 21-45 
110 | #440 | 0-057) 1,560 | 4,777 | 156 | 6,408 | 5,812 | 22-00 
120 | 6098 0-040) 1,760 | 6,556) 262 | 8,558 7,489 | 22-74 
1890 | 8-847 | 0-016) 1,060 | 8,981 | 11,810 9,508 | 23-76 
140 | 11-47 | 0-800/ 2,160 | 12,860 | 556 | 15,060 12,830 | 25-15 
150 | 15:89 |0-857/ 2,960 17,100 84% | 20,800 | 15,880 | 27-15 
160 | 22:89 | 0-819) 2,560 | 24,000 | 27,940 20,670 | 30-05 
170 | 82-40 | 0-764) 2,760 | 34,860 7 89,630 27,140 | 34-58 
180 | 40-19 | 0-686| 2,960 | 52,920 59,150 | 86,120 | 42-12 

26 


expressed as B.Th.U./(S.G.V.)(°F) and Aw as 
B.Th.U./lb. The slope is therefore (Ibs. 
H,0/S.G.V.)/(°F). Lines of constant wet bulb 
temperature also have the slope of — EA/y,. 
where E is the experimentally determined 
constant for the wet and dry bulb hygrometer. 

In the gas industry, the standard method of 
measuring the dew point of gas employs the 
G.L.C.C. hygrometer described in the Gas 
Examiners’ General Directions [11]. The 
instrument factor for this hygrometer at 30 in. 
Hg. and 82°F is given as 0-42, i.e. 


Pw — p = 0-42 (t — tw) (26) 


where p is the partial pressure of the water 
vapour in mm in this instance. Expressing the 
water vapour content as lbs. H,0/S.G.V. 


0-42 x 18 


a) = 762 x 0.885 (27) 


(Hw — (tg — tw) 


or Hy — Hg = 0-0258 (t, — tw) 


hg/kg Ay = 90-0258 
= 27-7 B.Th.U./(S.G.V.)(°F) (28) 


As aH is the correct potential for humidity, the 
wet bulb lines terminate on this saturation line. 
If it is desired to use values for the rates other 
than 27-7 B.Th.U./(S.G.V.)(°F), then constant 
transfer potential lines and wet bulb lines are 
no longer coincident. 


The enthalpy of moist gas may be expressed 
in the form 


ig = ¢@ (tg — to) + ey He (tg — fo) + Ao He (29) 


The three terms are the sensible heat of the dry 
gas, the sensible heat of the water vapour, and 
the latent heat of the water vapour, all related 
to the reference temperature (32°F). A 
constant value of 20 B.Th.U./(S.G.V.)(°F) is 
assumed for the specific heat of dry gas (cg) 
over the range considered, based on the com- 
position of Fulham mixed gas averaged over 
a year. The enthalpy at saturation 


is = ¢@ (ts — to) + ey He (ts — bo) + Hs 
In order to obtain a straight line plot on 
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the humidity chart, a term “ apparent en- 
thalpy " is introduced, defined as 


ig = Cgltg — to) + ev Hs (ty — to) + Hg (80) 
— ts = tg (tg — ts) + (Hg — Hs) (81) 
A line of “constant apparent enthalpy 
(ig — tig = 0) is thus a straight line of slope 
— terminating at the saturation line 
(Hg vs. ts) on the psychrometric chart. The- 


deviation between the true enthalpy and the 
apparent enthalpy is therefore given by 


ig — ig = ty Hg (tg — to) — ty Hs (ts — to) 


The form taken by equation (82) is indicated 
in Fig. 4. 


+20} 
20 80 100 
° tg -ts 
-20} 
Fic. 4. 


As ¢g and A, are both constants, all lines of 
constant apparent enthalpy have the same 


slope on the psychrometric chart, namely 
20 
— jor ‘= 00186) 


thus interpolation is a simple process. The 
difference between the true and apparent 
enthalpy is recorded for convenience on the 
chart as a series of curves of percentage 
deviation. It will be observed from Fig. 4 that 
it is possible to have three points where the 
deviation has the same absolute value, two 
being positive and one negative. 

Two psychrometric charts, illustrated by Figs. 
5 and 6, have been constructed in the manner 
described, covering different temperatures ranges. 
In calculating the enthalpy of moist gas a 
value of 0-45 B.Th.U./lb. for cy has been 
assumed throughout. This does not take into 
account the variations in specific heat with 
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HUMIO HEAT, S B.ThU/SGYV °F 
20-5 21-0 


yo 40 to 70 100 710 120° 
TEMPERATURE 
Fig, 5. Psychrometric chart for coal gas (at 30 in. Hg. 


pressure ), 


temperature and partial pressure, but these 
variations are within the accuracy of the 
charts. The most severe limitation of the charts 
is that they apply to a single total pressure, 
and other charts must be prepared for pressures 
other than 30 in. of mercury. 

The data used in the construction of the 


wUMIO HEAT, S 


MASS TRANSFER POTENTIAL, aH Ib water varour/SGV. 


charts and that required for the graphical 
method are summarized in Table 1. 


4. Aw ExampLe or tir Use or 
THE CHARTS 
Example 1-—Consider moist gas at 120°F with a dew-point 
of 60°F, entering a washer cooler. What is the maximum 
theoretical temperature to which the cooling water will 
rise ? 

The theoretical maximum temperature for the water is 
that at which the moist gas and the water are in equili- 
brium. Assuming that the ratio E for the washer-cooler 
is equal to that for the wet and dry bulb hygrometer, then 
the wet bulb temperature is the criterion. 

From the chart the wet bulb temperature for gas at 
120°F with a dew-point of 60°F is 83-0 F. If enthalpy 
were considered as the driving force then the limiting 
temperature would be that given by a constant enthalpy 
line, namely 78-8°F. 


5. APPLICATION OF THE GRAPHICAL 
Metuopv To CONDENSER DesiGn 


Po.LLarD [12] describes a test on a cross-flow 
water-tube condenser, and compares the area 
calculated by the method of and 
Hovcen [1] with the known area. Although 
the gas and water flown in such a condenser 
are not truly countercurrent, the example does 
provide an opportunity of comparing the 
proposed graphical method with the laborious 
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Fic. 6. Psychrometric chart for coal gas (at 30 in. Hg. pressure). 
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Example 2.— Cross Tube Flow Cooler. 


Gas Rate = 4,000 cu. ft. hr. at S.T.P. 
Cooling Range = 160°F to 53°F. 

Saturation Temperature of Inlet Gas = 156-5 °F. 

Water Rate = 2,000 Ibs. /hr. 

Water Temperature Range = 445° to 91-5° = 47 F. 

Enthalpy of Inlet Gas = 25,000 B.Th.U./S.G.V. 

Enthalpy of Outlet Gas = 1,100 B.Th.U./S.G.V. 
23,900 


= 458.G.V. hr. 


Slope of Operating Line 509 


B.Th.U.(S.G.V.)(°F). The values for the 
individual film coefficients used by PoLLarp 
are calculated values and allowance is made 
for varying conditions of gas and water flow. 
The proposed graphical method requires that 
mean values of the individual coefficients be 
used. The relationship between he, kg and h, 
and the nomenclature used by PoLLaRrD are 
as follows : 


Ig == hie, 


ke = and h, = hy 


2- 


The mean values for the individual coefficients 
based on those calculated by PoLLarp are : 


he = 10 B.Th.U /(ft®)(hr.)(°F). 
Ib. moles (atm. 
kg = 1-42 
= 0. 546 (ft®) (Ibs. (S.G.¥.)) 
h, = 80 (ft?) (hr.) (°F). 
E = he kg = 18:3 B.Th.U./(S.G.V.) (°F). 


and hy /kg = 147 B.Th.U.(S.G.V.) (°F). 


It will be observed that the value of E calculated 
by the j factor method of CurLton and CoLBuRN 
is approximately two thirds of the value obtained 
from the wet and dry bulb hygrometer. In order 
to obtain an exact comparison between the 
CoLsurx and Houcen method and the graphical 
method, an enthalpy transfer potential based on 
E = 18-3 B.Th.U. (S.G.V.)(°F) will be used. This 
function appears in Table 1, as column 12. 

It will also be assumed that the gas rapidly 
approaches saturation and remains in this con- 
dition throughout the condenser. This simplifies 
the problem considerably, as the enthalpy of the 
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gas and its enthalpy transfer potential are 
functions of the same saturation temperature and 
the use of the psychrometric chart is obviated, 
once the terminal conditions have been established. 
The validity of this assumption which is inherent 
in the and Hoverx method, will be 
discussed later. 
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Fig. 7 is the graphical solution of the problem 
with tie-lines drawn in at ten equal increments of 
gas enthalpy. The values of gas enthalpy, water 
temperature and transfer potential of the gas and 
the interface are given in Table 2. 

The value 1/0, — 0; for each tie-line is given 
in Table 3 and dig x 10% = 2.890 B.Th.U. - 


dig 2-39 
S.G. y Simpson’s Rule x 
10-61 = 8-45 Transfer units 
4 
A = 8-45 x —-— = 61-9 sq. ft. 


0-546 


G. 5. Criss and E, T. Newson 


Table 2. 


~ 10% 
B.Th.U./ 
8.G.V. 


(Og 4) 10% 
B.Th.U./ 
S.G.V. 


1 
(0g + Hg (A — As) 
x 
B.Th.U./ 
8.G.V. 


ag Hg (A — Ag) 
x 10% 
B.Th.U./ 
S.G.V. 


25-00 
22-61 
20-22 
17-83 
15-44 
13-05 
10-66 
8-27 
5-88 
3-49 
1-10 


18-90 
17-35 
15-90 
14°25 
12-50 
10-90 
9-07 
7-25 
5-29 
8-25 
1-05 | 


6-42 
5-29 
417 
3-12 
221 
1-43 
0-78 


01405 
O-1419 
O-1440 
01496 
01604 
0-1739 
0-1992 
02871 
0-3188 
0-5380 
3-651 


0-164 
0-166 
0-164 
0-164 
0-154 
0-140 
0-121 
0-097 
0-068 
0-037 
0-004 


6-95 
6-88 
6-78 
6-52 
6-08 
5-61 
4-90 
412 
3-08 
1-82 
0-27 


The value 1/@, — 6; for each tie-line is given in Table 8 and dig = 10° x 2.890 B.Th.U./S.G.V. 


Table 8. 


Value of the tie-line slope. 


— 147 B.Th.U./(S.G.V.\(°F) | 


«© — 80 B.Th.U./(S.G.V.\(°F) 


1/0g — 9; | 


0g — | 


0g % 


— 0; 


15-31 
14°22 
18-18 
11-85 
10-41 
9-10 
7:55 
5-06 
421 
2-36 
0-85 


6-95 
6-88 
6-78 
6-52 
6-08 
5-61 
400 
412 
3-08 
1:82 
0-27 


0-1453 
0-1497 
01534 
01645 
0-1783 
02041 
02427 
0-3247 
05495 
3-704 


0-2278 
02278 
02273 
0-2300 
02433 
02564 
0-2825 
0-3257 
04184 
0-6849 
4167 


00653 
0-0703 
00761 
0-0844 
0-0962 
0-1099 
0-1825 
0-1678 
0-2875 
0-4237 
2-857 


| 
1/0gq — 9% | 
| 


| 
| 
3:50 
O-414 


| 


5-90 
2-10 
0-856 


Ny = Number of enthalpy transfer units corresponding to total area. 
Ny = Number of enthalpy transfer units correpsonding to half the area (91-5°F to 58-9°F). 


This compares with 61-4 x 1-6 sq. ft. calculated 
by Pottarp. His addition of 1-6 sq. ft. for the 
range 160°F to 156-5°F is not valid because it 
assumes that no latent heat is transferred in this 
range. Referring to Fig. 7, the assumption is 
that the initial gas path is one of constant 
humidity (dH = 0). The initial slope of dig /dtg 


by graphical construction gives an enthalpy of 
28,200 B.Th.U./S.G.V. at 156-5°F, a decrease of 
1,800 B.Th.U./S.G.V. The amount of superheat 
in gas at 160°F saturated at 156-5°F is only 
100 B.Th.U./S.G.V. thus the correction should 

100 


be 


Ta00 * 1.6 or approximately 0-1 sq. ft. of 


tr ig x 10 x 104 
°F B.Th.U_/ B.Th.U./ 
S.G.V. S.G.V. | | 
91-5 1105 | 
86-8 10-47 
82-1 9-12 | 
17-4 7-73 | 
72:7 
| 
633 
58-6 
53-9 | 
49-2 
$45 | 
7" 
5 
 — 
915 | 440 
86-8 440 
82-1 | | | 440 
17-4 | 4-33 
72-7 #11 
68-0 8-90 
63-3 | 38-54 
58-6 38-07 
53-9 } 2-39 
49-2 1-46 
45 0-24 
| 
Nr | 10-88 
Na | oes 
Ny/Nr | | 0-464 
80 
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the 1-6 sq. ft. observed. The area therefore of 
61-9 sq. ft. calculated by the proposed graphical 
method agrees with that calculated by the 
method of CoLpurN and HovceEn (61.5 sq. ft.) to 
within 1%. 

It has been stated previously that the true 
driving force is given by the expression 


(9, — +(A — As) a5 Hs 


The magnitude of the deviation, the second 
term in the expression, appears at the sixth 
column in Table 2. The last column gives value 
of 1/[(0g — 4) + (Ay — As) Hy]. By Simpson's 
Rule 
(9g —9%) +(\ 8 
10.40 = 8-28. 


x 


4 

A = 8-28 x 0546 * 60.7 sq. ft. 
The effect therefore of omitting the deviation is 
to obtain an area 2% in excess of that obtained 
by using the true driving force. For design 
purposes the error is in the right direction, but 
when determining coefficients allowance should be 
made for the reduction in driving force. 

In order to check whether the ratio A, ‘ky 
assumed corresponds to the observed facts, a 
knowledge of the water temperature at the 
junction of the third and fourth passes (corres- 
ponding to half the transfer area) is required. 
Unfortunately PoLLarp does not give these data 
for the test considered, but from the profiles 
given in Fig. 1 of his paper, one may assume it is 
approximately 53-9°F. This temperature is 
selected for convenience in demonstrating the 
method, as it corresponds to eight out of the ten 
increments of gas enthalpy. Values of (@g — @,) 
have been obtained from Fig. 7 for two other 
tie-line slopes, and the number of transfer units 
corresponding to eight and ten increments 
calculated. The results are summarized in 
Table 3. 

The reasons for selecting the two other values 
for the tie-line slopes are as follows. Infinite 
slope implies zero liquid film coefficient and 
t, =t,. The number of transfer units may be 
calculated for this case without graphical con- 
struction of the tie lines. The latter value of 


— 80 B.Th.U./(S.G.V.)(°F) is assumed to be in 
the region of the correct value. 

The effect of plotting the ratio N,/N 7 against 
the reciprocal of the tie-line slope is shown in 
Fig. 8. The fact that the three points are colinear 
is fortuitous as this relationship cannot be 
established from equation (25). Assuming kg is 
a constant, the ratio N,,/N, should equal one 
half, if the tie-line slope is correct. If therefore 
58-9°F was the correct mid-point temperature the 
correct tie-line slope obtained by extrapolation 
would be — 60 B.Th.U./(S.G.V.)(°F). Normally 
the plot obtained will not be a straight line and 
the correct tie-line slope should be obtained by 
interpolation rather than extrapolation. 


Fic. 8. Correct slope of hy /kg by extrapolation = — 60 
B.Th.U./(S.G.V.)(°F). 


It is appreciated that kg cannot be a constant 
throughout the condenser due to the large change 
in humidity of the gas and thus in the true mass 
flow of the mixture. In this example the calcu- 
lated value of kg at the hot end is 1-28 times the 
value at S.T.P. If h, remains substantially 
constant, the effect of the increased value of k, 
is to decrease the tie-line slope and hence the 
value of (9g — @;). The total heat flux however 
depends on the product kg (@g — 8,) and the two 
factors compensate one another. The compensa- 
tion would be exact if the equilibrium line were 
of infinite slope. This is not so, but the slope 
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increases with humidity as does the mass flow, 
and the residual error is small. 

Fig. 9 is the graphical solution of the problem 
using an enthalpy transfer potential based on 
27-7 B.Th.U./(S.G.V.)\(°F) and the same values 
for kg and h,. The resulting number of transfer 
units is 7-70 corresponding to an area of 56-4 
sq. ft. The area is smaller, because the assumed 
transfer potential implies a heat transfer coefti- 
cient (hg) equal to 15-0 B,Th.U./ft.? hr. °F as 
against 10-0 B.Th.U. /ft® hr °F. 

It will be observed from Fig. 9 that the effect 
of taking E equal to 27-7 B.Th.U./(S.G.V.)(°F) is 
that the gas approaches saturation more rapidly. 
The construction indicates, in fact, that the gas 
becomes supersaturated, i.e. fog formation occurs. 
Under these conditions the basic differential 
equations no longer holds and all methods, in- 
cluding that of and HovucEen, become 
approximations. If the true value of E is as high 
as 27-7 B.Th.U./(S.G.V.)(°F), then the ratio b/« 
will exceed unity over the whole range considered. 
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Consequently the rate of mass transfer will be less 
than that required to keep the gas saturated, and 
fogging will occur with gas that is initially 
saturated. Water condensing in the bulk of the 
gas reduces the total quantity of heat to be 
transferred across the interface, and a smaller 
area is required. 

Even when the ratio 6 /< is less than unity, fog 
formation is still possible, because on a curve of 
rapidly increasing slope such as the saturation 
line, the slope of a chord is inevitably much less 
than the slope of the curve at the upper extremity 
of the chord. Insufficient importance has been 
attached to the formation of fog in cooler-con- 
densers, and care should be taken to see that 
coefficients determined under fogging conditions 
are not applied to the design of plant in which 
such conditions do not exist. 

Although the method has been described with 
specific reference to moist coal gas, it can be 
applied to any other system for which the 


requisite data are available or may be calculated. 
Furthermore the integral — 0, may be 
evaluated by the other methods described in the 


literature on simultaneous heat and mass transfer. 


Acknowledgements—-The authors wish to 
express their indebtedness to Mr. G. F. Epwarps 
for preparing the diagrams and assisting in the 
calculations and to the Gas Council and the 
North Thames Gas Board for permission to 
publish this paper. 


NOMENCLATURE 

A = area of transfer surface 
E = ratio of transfer coefficients hg/kg 
(S.G.V./hr] 
[Ibs /S.G.V.) 


(ft?) 


G = gas flow 
H = humidity 
K = overall coefficient of enthalpy transfer 
(Ibs. 
= mass flow of liquid [Ibs. /hr} 
N = number of gas film enthalpy transfer units 
P = total pressure {in. of mercury} 
b = Lewis Number hg/kgs 
= heat capacity [B.Th.U./(S.G.V.)°F)) or 
(B.Th.U. (ft? 
{B.Th.U./S.G.V.} 
(B.Th.U./S.G.V.) 


h = heat transfer coefficient 
i = enthalpy 
i = apparent enthalpy 
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k = mass transfer coefficient (Ibs. A = latent heat [B.Th.U. /Ib.} 
or enthalpy transfer coefficient [B.Th.U./(ft®)(hr)( 
= slope of saturation lined@, 
= partial pressure of diffusing vapour (in. of mercury] G = gas or gas film coefficient 
= heat flux ([B.Th.U. /(hr\(ft®)} = interface 
humid heat cg + 0-45 11 (B.Th.U. (S.G.V.)CF)} LL = liquid or liquid film coefficient 
= temperature (°F) S = saturation 
humidity correction factor V = vapour phase 
fractional distance thourgh the gas film W = wet bulb conditions 
enthalpy transfer potential [B.Th.U./(S.G.V.)( F)] O = reference conditions (32°F) 
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Asbtract—The separation of common gases in pilot plant thermal diffusion columns has been 
investigated. The effects of the variables temperature, concentration, flow conditions, and 
column dimensions have been obtained experimentally. 

Results show that the equations deduced for isotopes apply to common gases in certain limited 
cases : that turbulence reduces enrichment, but that it is offset to some extent by a probable 
effective reduction in annulus width. Continuous operation of columns demonstrates a marked 
adverse effect of increase in withdrawal rates, which have to be greater than those normally 
used for isotope separation, if the method is to have industrial importance. Consequently, the 
column separation performance factors H and K, normally assumed independent of the small 
withdrawal rates for the case of isotopes, are no longer so for common gases. It is apparent 
that, for the separation of common gases, pilot plant design data is essential, and that a semi- 
empirical approach is necessary in the absence of theory relating the withdrawal effects. 


Résumé—Les auteurs ont étudié la s¢paration des gaz usuels dans une installation pilote de 
colonnes de diffusion. Ils ont étudi¢ expérimentalement l'influence des variables : température, 
concentration, conditions d'écoulement, dimensions de la colonne. 

Les résultats montrent que les équations valables pour les isotopes s'appliquent aux gaz 
usuels dans certains cas limites: a savoir que la turbulence réduit l'enrichissement, effet 
compensé dans une certaine mesure par une réduction efficace probable de la largeur annulaire. 

Les opérations continues mettent én évidence un effet notable antagoniste de l'augmentation 
de la vitesse de soutirage : pour que la méthode ait une importance industrielle, cette vitesse 
de soutirage doit étre supérieure aux vitesses normalement utilisées pour la séparation des 
isotopes. Par suite, les facteurs H et K d'efficacité de fonctionnement des colonnes qu'on suppose 
habituellement indépendants des petites vitesses de soutirage dans les cas des isotopes ne le sont 
plus dans les cas des gaz usuels. 

Il est évident que, pour la séparation des gaz usuels, il est nécessaire de disposer de données 
d installations pilotes et de procéder semi-empiriquement en l'absence de théorie relative aux effets 
de soutirage. 
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INTRODUCTION 


Wuews a gas mixture is passed through an annulus 
across which there is a temperature difference, a 
separation of the components of the mixture 
occurs by thermal diffusion. The effect was anti- 
cipated by Cuapman [1] and Enskoo [2], and 
considerable experimental work by Isrs [3] 
substantiated the theory. 

The effect is considerably enhanced by the 
convective flow of the gases on the walls of a 
vertical separation column. This was demon- 
strated experimentally by Ciustus and Dicker 
[4], who also proposed a preliminary analysis of 
the hydrodynamics of such columns. These 
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columns have received considerable attention for 
the purpose of isotope separation. 

The theoretical hydrodynamics of the process 
has been examined by CLusius and Dicket [4), 
Jensen [5], Stetnweve [6], and Watpman [7] 
of the German school, and by Furry [8, 9}, 
Onsacer [8], and Jones [8, 9] of the American 
school. The principal references mentioned 
contain a large number of additional references 
to isotope separation. The important features of 
the design of columns for such separations are the 
similarity of the molecules concerned, and the 
exceedingly small withdrawal rates used. 

There is little information on the application 
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of the existing theories. which were specifically 
derived for isotopes, to the design of columns for 
the separation of industrial gases. 

The following results of a pilot plant investiga- 
tion of the separation of mixtures of hydrogen, 
carbon monoxide, carbon dioxide, and nitrogen 
by thermal diffusion, give an indication of the 
complexity of the design problem. 


GENERAL PRINCIPLES 


For total reflux or “ discontinuous operation. 
the net differential transport of the heavy com- 
ponent in a binary mixture is given by the 
differential equation 


r=H-c(l —c)—K-de dz (1) 


The re-mixing effects opposing the separation 
by thermal diffusion are represented by the term 
AK - de dz. Furry further defined the individual 


re-mixing effects. 


(i) convective re-mixing, k, 
(ii) re-mixing due to back diffusion along the 
column, kg 
(iii) re-mixing due to asymmetrical effects. k,, 


Then 
K =k, +k, 


Under total reflux conditions the net differentia! 
transport of the heavy component r, is zero. Hence 
between appropriate limits, 


L c, (1 — co) 
N= gm ~ —<,) 
where K H = $A. and corresponds to the “ ideal 
separation length’’ proposed by the German 
school. Then N becomes the number of ideal 
separation lengths corresponding to a column 
height L, under given conditions of concentration 
and temperature. In the following work, the 
results are calculated with common logarithms. 
so that N’ = N/2-808. 

According to the theory, H and K are constants 
for any apparatus, their values being determined 
by the properties of the gas mixtures, and the 
dimensions of the separating columns. For 
simplified conditions, which apply to the annular 
columns in the following experimental work, 
Furry proposed the use of the “ plane parallel 


case. The properties of the gas mixtures are 
taken at the arithmetic mean temperature of the 
annulus. The expressions for the constants 
proposed by Furry are then as follows. 
2 
96n T 
k _ B gow Fl | 
2304-7? - D 
= 2-w-B-p-D > 


These expressions are developed on the assumption 
that no convective turbulence exists in the annuli. 
For the truly plane case, Furry suggests that 
turbulence occurs when Re > 25. and for safe 
working Re < 10. 


AT 


It will be evident from the equations that 
transport is aided by large values of w and B, 
but at the expense of increase in re-mixing effects, 
and decreases in separation. This is characteristic 
of thermal diffusion columns, necessitating a 
compromise which can be made eventually only 
on a basis of experimental work. 

For continuous operation, of a column where a 
definite amount of concentrated material is with- 
drawn continuously, the differential equations of 
STEINWEDEL or Furry may be used. The 
integrated STEINWEDF!. equation has been pre- 
ferred for this work. 


where 


y= oH and T- 


To design a thermal diffusion column for 
separating a certain mixture, it would be necessary 
to establish : 


(i) that existing hydrodynamic equations 
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satisfactorily correlate experimental re- 
sults. 

(ii) the values of column transport coefficients 
for different sizes of columns under differ- 
ent conditions, and to show the effects of 
these variables on H and K. 


(iii) to what extent conditions approach turbu- 
lence without being too harmful to the 
separation. 

(iv) how the constants or coefficients H and K, 
once determined, can be used for scaling 


up plant. 


Propos—eED METHOD FOR DETERMINING THE 
EXPERIMENTAL VALUES OF THE COLUMN 
Transport Factor H UNDER GIVEN 
CONDITIONS 


A theoretical plot of the total reflux conditions in 
a column can be represented as in Fig. 1 (eq. 2). 
For the pilot plant column of length L, the ideal 
separation length }4 or K/H would be fixed by 
the physical conditions and the dimensions of the 


apparatus. Accordingly, N and N’ will have 
practical values for the set conditions. The value 
of N’ can be established by plotting c, versus cy 
from the experimental results of total reflux 
operation. As L for the column is fixed, N’ yields 
an experimental value for K/H. If now an 
experimental value for H could be obtained for 
the same column, then an experimenta! value 
of K could be deduced. 

The STEINWEDEL equation for continuous 
operation can be represented graphically as in 
Fig. 2 with constant initial concentration of heavy 
component ¢y. For a practical column with fixed 
initial concentration co, the end-concentration c, 


corresponding to a series of withdrawal rates o, 
can be determined experimentally. Under con- 
stant temperature conditions, a curve similar to 
that shown in Fig. 8 could be obtained from the 
results. 

The K/H and N’ values for the column are 
available from total reflux experiments. Choosing 
a value of o within the experimental range, say 
o,, the experimental end concentration corres- 
ponding to cg would be c, from Fig. 3. From 
Fig. 2, N’ and cy being fixed, the value of y could 
be determined which gives the correct ¢,, viz. ¢,,. 
Suppose the value of y is y,, then y, = o,/H or 


Fic. 8. Constant T, AT, C, and L’. 


H =0,/y,. In this way the value of H corres- 
ponding to the experimental condition is deter- 
mined. 

Having obtained (K/H),,,, and H,,,,. and by 
deduction K,,,,, their values can be compared 
with theoretically calculated values. Such a 
comparison should indicate to what extent 
design could be fundamental for the case of 
common gases. 


DESCRIPTION ON THE PILOT PLANTS 


Two plants were built. The first consisted of a single 
column 15 ft. high, gas heated, and the second of three 
6 fi. electrically heated columns. The columns differed 
in their diameters and annuli widths. The dimensions 
are given in Table 1. Line diagrams of the columns, in 
Figs. 4 and 5, describe the construction. High water 
rates were maintained on the outside of the outer pipe 
for cooling, and the upper distribution ring was swathed 
in muslin to reduce splashing. 

Needle valves on the separate heating-gas lines of the 
15 ft. column permitted good control, and a reasonably 
uniform temperature along the length of the column. 

The electrical heating of the 6 ft. columns was con- 
trolled through external resisters. The heaters were of 
nichrome wire threaded with recessed porcelain beads. 
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Connections were at the top only, and the wire and beads 
were plaited before inserting into the central tube. Tem- 
perature measurements on the periphery of the inner 
tube were not taken, so that it is not known to what 


= 
MILO STEEL 


gORAIN 


AIR INLET FOR BURNERS 
Fic. 4. 15 ft. column. 


extent temperature assymmetry interfered with the 
separation. The effect would not be important for the 
two larger 6 ft. columns, The inner tubes of the 6 ft. 
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Fra. 5. 6 ft. column: aa—Circulating gas connections ; 
bb —Product gas connections, 


columns were of stainless steel, and it is now realized 
that for the smallest column the preferred material would 
be copper, as, with its much higher thermal conductivity, 
asymmetrical temperature effects wou!d be minimized. 


,Diametrical asymmetry was minimized by fitting a 


centralizing spider to the outside of the inner tube. The 
material used in the construction of the top flanges and 
the outer tubes was heavy gauge mild steel to obviate 
buckling under heat stress. 


I 
ils 


GAS - HOLDER 


Fic. 6. Thermal diffusion pilot plant : N - needle valves 
S - sample cocks ; Tr.C three-way cocks. 


A line diagram of one of the pilot plants is shown in 
Fig. 6. The columns could be operated in parallel or 
singly, and a re-arrangement was made later to operate 
in multistage. 

Photographs (7 und 8) showing some detail and a gen- 
eral view of the plant are given in Figs. 7 and 8 respectively. 

“ Rotameters " were used to measure gas rate. They 
were calibrated for each mixture over the concentration 
range being measured, by means of a soapfilm bubble 
flow-meter. 

Gas analysis was by means «f a katharometer, and by 
direct Orsat. Small samples were taken, or the sampling 
spread over a period to cuuse least disturbance of the 
conditions. 

For all cases, gas was circulated around the tops of 
the column at a relatively high rate (approx. 20 1/min.). 

Continuous operation necessitated a balanced operation 
as the system was sensitive to small volume and pressure 


changes. 


EXPERIMENTAL RESULTS 


These were obtained for the following mixtures. 
H,/CO,, H,/CO, CO/CO,, and O,/N, in the three 
6 ft. columns, and for H,/CO mixtures in the 
15 ft. column. 

The results cover both the “ total reflux” or 
“discontinuous operation,” and continuous 
operation " having a finite fixed withdrawal rate. 


COOLING 
CIRCULATING | 
bb Proouct | ond 
— 
| 
| 
87 
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EXP VALUES 


41=250°C 


020 40 80 10 
Co “le 
Fic. 9. Total reflux Hy-CO. 15 ft. column. Fic. 10, Total reflux H,-CO,. 


Column 2. 


1000 
CCS/MIN 


Fic. 11, Continuous operation Hg-CO,. Column 2. 


Column dimensions. 


Co*220%%, CO, 
1} © wO 
x 
CCS/MIN 


Fig, 12. 


Continuous operation Hy-CO,. Column 1. Inner tube o/din. | 
Outer tube i/d in. | 4° | 8° 


poe 


em.) 2805 | 1910 


| 


ft. 


6 | 6 
220% CO, 
41 #280 °C 


. | These are recorded in Tables 2—9 inclusive. 
_— Certain of these results are plotted in Figs. 9—13 
O 160° 260 300 400 S00 inclusive. The corresponding tables and figures 
are as follows :—Table 2(b), Fig. 9; Table 2(a), 
Fig. 10; Table 5. Figs. 11, 12 and 13. 


Fic, 14. Continuous operation Hy-CO,. Column 3, 


pnd 30 THEOR CuavE 
| 
x~ EXP VALUES 4 
at 20° | 
2 ? 
10 
2 
16 60 
vol. 
1956 
6 
d Tale 
‘ 
| Columns | 
ac’ 
| 64 ‘32 ‘32 
1-80 | O48 | 215 
— 
| 8" /gs 
[ 1460 | 42-80 
6 15 
30 
20 
|. 
38 


Fic. 7. A general view on the 6 ft. thermal Fic. 8. A detail view around the main 
diffusion columns. blower of the thermal diffusion plant. 


} 
OL 
| | | | | i 


1956 


VOL. 
5 


The separation of common gases by thermal diffusion 


Table 2 (a). Discontinuous or total reflux operation H,,CO, mixtures. ¢ — \ CO, by volume in the mixture. 


| 


T,°C AT | [K/H less. 


250 
800 


100 
155 
200 
250 
250 
250 
250 
250 
250 
250 
250 


250 


Table 2 (b). % CO by volume in the mizture. 


TC [K/H expt. 


270 735 
320 516 


120 


1 | 3-58 22-0 26-4 270 | | | 7174 0-105 
| 3-58 22-0 27-5 320 | 615 0-129 
2 1-80 22-0 29-0 Ca 60 | 505 | O-157 
1-80 22-0 31-8 120 368 | 0-216 
1-80 22-0 | 35-5 | 175 276 | 0-288 
1-80 22-0 227 | 850 
1-80 22-0 41-2 | 270 205 0-388 
1-80 5-2 12-0 | 270 200 0-400 
1-80 10-0 23-2 | 270 183 0-434 
Fig. 10 1-88 12-5 27-0 270 192 0-412 
1-80 16-0 82-1 | 270 202 0-394 
1-80 22-0 40-8 270 190 0-388 
1-80 25-2 45-1 | 270 190 0-388 
1-80 85-2 56-0 270 216 | 867 
| | 
5 3 0-48 22-0 63-0 | 270 = 102 0-790 
956 
| | | | 
Ne | | | 
1 3-58 15-0 18-6 0-110 
3-58 15-0 | 20-1 | 0-154 
2 1-80 15-0 22-7 = 100 | 360 0-821 
1-80 15-0 | 220 20000 | 
1-80 15-0 20-1 270 250 217 0-366 
1-80 15-0 30-8 820 300 198 0-404 
3 0-48 15-0 30-6 | 270 250 | 200 | 0-398 
0,48 15-0 31-6 320 300 190 0-418 
0-48 15-0 82:5 875 355 184 | 0-434 
15 ft. 2-15 50 35-0 215 200 22 | 1 
215 6-2 43-3 215 | 200 191 1-06 
2-15 0-8 53-8 215 | 200 193 | 1-05 
Fig. 9 2-15 10-4 54-0 215 200 192 1-058 
2-15 13-5 60-5 215 200 194 1046 
215 21-0 67-0 215 200 230 | (0-882 
Table 3. Discontinuous operation. H,/CO (water-gas) 15 ft. column. 
CO, 2-60 3-30 
co 47-30 61-8 66-30 89-6 
Ne 11-90 20-00 
H, 38-20 10-40 
39 
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Table 4(a). Discontinuous operation. O, Ng. 


| 


1:80 210 224 270 250 2310 0-035 
| | 


j | | 

| O48 21-0) 266 320 300 598 185 

j 
| 


| 


Table 4(b). Discontinuous operation. CO,CO,. 


Col. 


2 180 320 33:5 270 250 2700 


048 32:0 462 270 250 
| 


Table 5. Continuous operation Hy CO, ¢ = °, by volume. Withdrawal rate = a ce min. 


T, = 20°C, cg = 22-00%, COg T, = 270°C, Column 2, Fig. 11. 


38-4 36-1 34-2 43-6 32-0 
220 580 640 


T, = 20°C, cg = 22-00% COg, T, = 220 C, Column 2, Fig. 11 
220 


T, = 20° . 175°C, Column 2, Fig. 11 


310 
380 


T, = 20°C, cg = 220% T, 
28-4 274 26-3 25-2 
310 520 800 1280 


T, = 20°C, cg = 220% COg, Ty = 80°C, Column 2, Fig. 11 


26-8 25-6 248 24:3 
220 480 820 1120 


T, = WC, cg = 22-00% COg, T, = 270°C, Column 1, Fig. 12 
25-9 25-6 25-3 25-2 25-2 
220 350 510 730 1120 


T, = 20°C, cg = 22-00% CO,, T, = 320°C, Column 1, Fig. 12 
25 26.3 26.1 260 
330 450) 620 Wo 


T= 20°C, = 220% T, = 270°C, Column 3, Fig. 13 
33-5 26-0 242 24-0 23-6 
110 170 230 250 400 


40 


Col. 
(0090 
| 08 0-262 
| VOL, 
5 
195¢ 
41-2 20-7 
30-0 28-6 
2 | SS... = 
13-6 29-4 28-0 26-4 
2, Fig. 11 
31-8 
200 
264 
27-5 
63-5 
0-0 


Table 6. 
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Continuous operation H, CO. 


c 


% by volume. 


T, = 20°C, T, = 320°C, ¢, = 15-0%, CO, Column 2 


27-7 


27-0 


25-7 
660 


860 


23°3 
1120 


270°C, 


23-7 
630 


220°C, 


21-2 


= 15-0% CO, 
22-9 
786 


Column 
22-4 
880 


15-0°% CO, 
20-6 
880 


Column 2 
28:2 
1000 


20-4 
400 


= 170°C, 


19-2 
620 


810 


15-0% CO, 
17-7 
1120 


Column 
17°5 
1260 


174 
200 


170 
270 


T, = 20°C, Ty = 270°C, 


16-9 
500 


15-0% CO, 
16-7 
1020 


Column 1 


T, = 20°C, T, = 220°C, ¢ 


17-9 


17-7 
650 


1200 


15-0% CO, 


Column 1 


T, = 20°C, T, = 270°C, cy = 


15-7 


15-7 
100 


15-7 
200 


CO, 


Column 3 


17-2 
45 


17-6 
100 


T, = 20°C, Ty = 320°C, 


16-8 
135 


15-0% CO, 


Column 3 


T, = 20°C, Ty = 875°C, ¢ = 


17-6 


| 


17-2 
175 


15-0% CO, 


Table 7. 


Continuous operation O,/N, 


21-46 
200 


Column 2 : 


T, = 20°C, T, = 270°C, ey = 21-00% O, 


21-1 
179 


T, = 20°C, T, = 820°C, cy = 21-:0% O, 


Continuous Operation CO/COs 


Column 2: 


T, = 20°C, T, = 270°C, cy = 320% CO, 


| 


33-0 
310 


33-0 
4+) 


82-6 
200 


T, = 20°C, T, = 270°C 


» = 820% CO, 


41 


30-8 28-7 | 22-1 21-3 
| 360 | o | 1420 | 1700 | 
T, = 20°C, T, = 
20-1 26-4 25-3 24-6 21-9 21-1 20-0 
o 00 | 240 360 480 1000 1230 1649 
T, = 20°C, T, = 
27-0 24-7 23-7 22-5 19-2 18-8 18-7 
e 0-0 200 330 500 - 1420 1700 1800 
T, = 20°C, T, - 
22-7 21-3 17-3 
_ | | nana | | 
18-5 17-9 
20-1 18-6 
30-6 15-8 
318 19-6 
be 
Colum 
ts | 32-7 | 18-3 | | 
00 50 
G | 22-4 | | 21-5 | 21-4 | 21-4 | 
0-0 250 820 580 
Column 
G | 26-6 | 21-2 | 21-1 | | 
e 0-0 (0 100 
G 33-5 33-1 33-0 33-1 82-9 42-8 32-8 
e | 0-0 | 220 | 220 | 260 | 380 Se 640 800 
Column 
% | 46-2 | 83-8 | 83-0 | | 82-6 | 
e 0-0 75 105 400 
|_| 
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Table 8. Continuous operation H, CO. 


15 ft. column. 


T.= 20°C. T, = 270°C. 


f = 50% CO 


15-7 
800 


16-5 
2360 


Co = co 


18-3 
1400 


1500 


= 97°,CO 


25-8 
1060 


28-1 
2180 


= 141% CO 


346 
800 


Table 9. Power consumption. 


6 ft. columns. 


H,/CO, : 


Column 1 


f = 220% COg; T,= 20°C; Ty = 270°C 
5 amps, 155 volts across the heater, i.e. 775 watts. 


On the outer pipe surface : 0-12 watt/sq. cm. or 0-775 watt./sq. in. 


Column 2 


6-2 amps, 82 volts, i.c. 500 watts. 0-10 watt/sq. cm. or 0-635 watt/sq. in. 


Column 3 


620 watts, 0-17 watt/sq. cm. or 1-10 watt/sq. in. 


Analysis of experimental results 

The results are re-expressed below in suitable 
units for intrepretation. The values of D, the 
diffusivity for various mixtures at different 
temperatures were obtained from the semi- 
empirical equation of GiLLILanp [10). Densities 
of mixtures were obtained at appropriate mean 
temperatures, and the arithmetic mean of the 
concentrations at the ends of the columns were 
used to establish a mean density. For the sake 
of simplicity such a procedure was considered 
adequate at this stage. The viscosities of gas 
mixtures were taken as being proportional to 
concentration for the same reason. Over the 
concentration ranges being considered the errors 
introduced should be small. 


The values of KH from experimental data. and a 
consideration of turbulence in the columns 


From the results of total reflux operation of the 
columns, it would be useful to determine how the 
‘ideal separation length’ (on the Jensen basis), 
or the ratio K/H (on the Furry basis), is affected 


by temperature and concentration for any 
mixture. The theoretical values of KH and Re, 
for the experimental conditions have been 
evaluated from the following equations, and are 
given in Table 10 for the separation of a H, CO, 
mixture in Column 2. 


K 084 T [Ar* 


AT 


(JENSEN) 


In order to obtain the values for [A /H ]tneor.- it 
is necessary to choose a value for x. This is 
extremely difficult as WALDMAN [11] gives values 
varying from 0-17 to 0-41, with concentration 
changes from 0-96 to 0-081 mol. fraction of CO, 
in H,/CO, mixtures at 20°C. A mean value of « 
of 0-25 was taken as being reasonable over the 
concentration range required. WinTeER [12] gives 
values of x at two temperatures, but while it 


G 90 13-1 16-5 35-0 
Fo 13-4 61 20-0 22:5 
5-9 17-7 21-7 26-4 33-3 39-4 
188 22:7 9 41-6 41-9 610 
VOL, 
1956 
3 
AT 
u*.g-p* AT 
42 


956 
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Table 10. 


Total reflux or discontinuous operation H,/CO,. cy = 22-0% CO,, Column 2. 


Ce expt. AT 


[K/H expt. [K/H]theor. 


368 
276 
227 


| 
| 
| 


demonstrated that a is temperature dependent. 
no correlation is possible. Nrer [18] obtained a 
relationship « « T~* for Neon isotopes at low 
temperatures. On the basis of present information, 
it is unsafe to take such a relationship. A plot 
of K/H versus T /AT gives similar curves for the 
experimental and theoretical case. As to the 
condition of the convective flow, it would appear 
that in every case except Column 3, turbulence 
exists, i.e. Re, > 25. This is shown in Table 11. 


Table 11. Total reflux or discontinuous operation. 
,/CO,. = 20°C. c= % co,. 


expt. 


80 
120 
195 
220 
270 
270 
270 


It will be observed that, despite the fact that 
convective conditions are turbulent, the experi- 
mental values all lie on the appropriate equilibrium 
curves. This suggests that the separation is not 
interfered with to the extent which would be 
expected due to re-mixing. This applies to Re, 
values as large as 160. An explanation would be 
that the temperature gradient is reduced in a 


central region between the hot and cold walls, 
due to a higher heat transfer rate in the central 
turbulent core, so that the existing gradient 
occurs across a reduced width. A smaller effective 
annulus width would result in an increased 
separation, which would then tend to counteract 
loss of separation due to the re-mixing effects of 
turbulence. 

The results for the total reflux operation with 
H,/CO mixtures are given in Table 12. 

For the 15 ft. column the N’ value is reasonably 
constant at constant temperature over a wide 
range of concentrations. An interesting feature 
of the calculations is that, for columns 2 and 3, flow 
conditions are substantially streamline. It 
would then be expected that [K/H].. would 
deviate less from [K/H ]tneor, and it should be 
useful to obtain values of x from the JENSEN 
expression for K/H, but using [K/H],,,». values. 
Table 12 summarizes the flow conditions and the 
a values obtained. 

Watpmawn [11] has collected the results of 
Inpps, Grew and Hirst [14] for H,/N,. As 
nitrogen is almost identical with CO in its thermal 
diffusion effect, these values are of interest for 
comparison. They are given in Table 13. 

yN, = 0-20 and « = 0-33 from these collected 
values, whereas from Column 2 (Table 12), 
a = 0-148. If the WaLpMAn value is correct an 
assumption of « = 0-33 for substitution in the 
JENSEN equation to calculate [K/H] would give 
[K/H) values different from those obtained 
experimentally. The « values obtained for the 
15 ft. column are of the same order as those 
quoted by WaLpMan, despite the fact that 
turbulent conditions exist. 


| | Receate 
29-0 80 60 | 570 | 47-4 
31-8 120 100 492 66-5 
35-6 195 175 386 82-5 
39-0 220 200 ) 329 83-5 
41-2 270 250 205 | 287 86-0 

OL. 
5 

Col. No. Re, 

2 22-0 29-0 47-4 

81:8 66:5 

- 85-6 82-8 

39-0 83-5 

- 41-2 86-0 

2 22-0 420 88-6 

0-0 60-0 - 160-0 

5-0 12-0 - 21-0 

43 
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Table 12. 


Column 2, Hg/CO ; eg = 15% CO; T, = 20°C 
220 270 | 320 


Total reflux or discontinuous operation. 


Column 3: H,/CO; cg = 15% CO; T, = 20°C 


15 ft. Column : H,/CO; T, = 15°C; Ty = 215°C 
6-2 9-3 10-4 
53:8 54-0 
188 197 
61:8 67-6 
1-056 1-005 
0-406 0-393 


43-3 


0-360 | 


Table 13. Collected values for H,/N, miztures. 


0-05 0-075 0-15 
0-47 0-42 0-40 


Y Ne 
0-47 


220°C 


0-75 
0-23 


0°55 
0-27 


0-85 
0-32 


0-25 
0-32 


The conclusions that can be drawn appear to 

be as follows, 

(i) the JENseN type of equation for [K/H] 
applies. 

(ii) It is impossible to predict with certainty 
from theoretical equations, the value of 
[K/H] as « is concentration and temper- 
ature — dependent, and the turbulence has 
an effect which cannot be correlated by 
existing equations. 


It is reassuring to find that, with fixed 
temperature conditions as for the 15 ft. 
column (Table 12), the [K/H] value 
remains reasonably constant over a wide 
range of Re, values, although it is under- 
stood that the factors influencing enrich- 
ment must be changing. Hence, it can be 
assumed that the only values of [K/H] of 
use for design are those determined experi- 
mentally, and which are relative to the 
conditions of the experiment. 


Tue or From 
EXPERIMENTAL DatTa 


The induction method used for obtaining 
experimental values of H has been described. 
The physica] data required for calculation of the 
values were evaluated at the appropriate arith- 
metic mean concentration relative to each 
experimental run. The conditions and experi- 
mental values of H are given in Table 14. 

It should be noted that, in the tables, o the 
withdrawal rate is given in cc/min as directly 
measured. For the calculation of H from experi- 
mental data, H = o'/y, the mass flow o’ g/sec, 
being calculated from the volume flow under the 
appropriate conditions. 


Tue Vatues or K OBTAINED DIRECTLY 
FROM THE EXPERIMENTA! DATA 


From total reflux or discontinuous operation, 
values of [K/H],x, have been obtained for the 


columns. 


Ce expt. 22:7 27-0 29:1 3-08 
ape. 360 247 217 198 
| 0-322 0366 0-401 
x | O-114 0-107 
T,°C 270 320 | 
30-6 31-8 
200 19 
N’ 0308 | 0-420 
Re, 046 - 
a 0358 0-360 
Ce expt. 35-0 60:5 67-0 
[K/H 195 202 
Re, 39-2 79-0 96-0 1956 
N’ 101 0-981 0-883 
0-464 0-412 
7 0-85 
0-21 
: 
44 


956 


Table 14. 


Column 2 cy = 22-0% CO, « = 0-25. 


T = 50°C, AT = 60, [K/H),,.4. = 505, Re, = 474 


27-4 
0-0024 
150 


28-6 
0-0024 


T = 70°C, AT = 100, xp. = 368, Re, = 66:5 


28-85 
0-C020 
250 


30-60 
0-0011 
70 


- 


T = 97-5°C, AT = 155, [K/H] 454. = 274, Re, = 828 


28-85 
0-0035 
670 


30-90 
0-0030 
890 


| 33-10 
0-0025 
160 


85-18 
0-0016 
20 


°C, AT = 200, [K/H),.5¢. = 227, Re, = 83-5 


29:8 
0-0046 
920 


82-4 
0-0089 
520 


35-50 
0-0037 
250 


88-00 
0-0045 
60 


°C, AT = 250, [K/H] exp. = 205, Re, = 86-0 


29-9 
0-0052 
1220 


83-5 
0-0042 
640 


86-75 
00050 
850 


| 
oe 


The values of H,,,., are obtained as shown for 


the same experimental conditions. 
values are directly obtainable. 
below in Table 15. 


Table 15. 


These 


Hence 


are given 


Column 2, ¢y = 22-0% CO,, = 0-25. 


T = 50°C, AT = 60, [K/H),4, = 505, Re, = 474 


K 


| 26-9 
0-57 
210 


27-4 
0-61 
150 


| 28-6 
0-326 


= 70°C, AT 


= 100, [K/H), = 368, Re, = 66-5 


26-9 
1-225 
630 


28-85 
0-728 
250 


0-408 
| 70 


30-60 | 


T = 97:5°C, AT =1 


55, [K/H] exp, = 276, Re, = 82-8 


Kept. 


28-85 
0-960 
670 


80-90 
0-836 
890 


83-10 | 
0-687 
160 


35-18 
0-427 
20 


= 120°C, AT = 200 [K/H] 


¢. = 227, Re, = 88-5 


29-80 
1-040 
920 


82:40 
0-878 
520 


85-50 
0-845 
250 


88-00 
1-015 
60 


°C, AT = 250, [K/H] x5, = 205, Re, = 86-0 


29-90 
1-055 
1220 


83-50 
0-865 
640 


86-75 
0-995 
850 


89-80 
1-420 
100 


The separation of common gases by thermal diffusion 


THEORETICAL VaLues or H 


These values have been calculated from the 
Furry equation for the experimental conditions 
in Column 2 with H,/CO, and H,/CO mixtures 


as shown below in Tables 16 and 17. 


‘gw 
96 


(Furry) 

‘The physical properties of the mixtures such 
as p and », were evaluated at the arithmetic mean 
concentration over the column for experimental 
conditions. 

As H is proportional to a, any change in « due 
to concentration or temperature will be reflected 
in H. From Wa.pMan, it would seem that 
x = 0-25 is a good average figure to use over the 
concentration range of the experimental case for 
H,/CO, mixture. The values of Hy... would 
then have to be used in the knowledge of its 
variation with concentration and temperature. 
The former can be related to « reasonably well, 
but the temperature effect is not certain. 


Table 16. 


Column 2, = 22.0% = 0-25. 


om 


T,°C 
AT 
H heor. 10* 


Similar values of H are calculated for Column 2 
on H,/CO mixtures on the basis of « = 0.33 from 
Wa.pman’s collected figures. These are recorded 
in Table 17. 


Table 17. 


Column 2, H,/CO, = 15-0% CO, = 0-33. 


| 18-85 


210 | 
26-9 
| 
| 
H expr. | | | 
- | | | 
- 
| 
 - | 80 | 120 | 175 | 220 | 270 
MB apt. | | | 60 | 100 | 155 | 200 | 250 
24 | 54 | 184 | 185 | 174 
| 
| 
— 
| 
|= 
| Cm 210 | 22.05 | 22-90 
120 22000 270 | 
| | AT 200 250 
H theor, X | 198 | 498 | 568 | 71-4 
45 
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Table 18. Total reflux operation. 


Column 2: Hy/CO,; cy = CO,. 


20-0 81-8 346 

TC 120 

AT 100 

ke 3-65 

kg» 105 

Re, 66-5 
Ree 9-56 


[ke /ka}* 
A 


| 
| | 


60-0 


41-2 
| 270 


89-0 
220 


| | | | 


THRORETICAL K 
These have been evaluated from the Furry 
equations below, using the experimental data on 
Column 2 for H, CO, and H, ‘CO mixtures, under 
total reflux. 


VaLuEs oF k, kg 


(Furry) 


9804 


Kk, + kg (neglecting k,) 


The values are given in Tables 18 and 19 below. 


Table 19. Total reflux operation : 
Column 2 : H,/CO; = 15.0%, CO 

30-8 
820 
300 
0-593 
12 
23-5 
3-27 


22:7 27-0 


220 
200 250 
0554, 0-592. 
100 | 107 
191 
0-564 0-593 


| | 


It is seen that k, the re-mixing effect due to 
convection is much greater than k, the longitu- 
dinal back-diffusion. For isotopes, Furry showed 
that in the streamline range Re,/[k,/kg}* is 
reasonably constant at 2.0. Both K, and kg 
increase with Re,. As given in Table 18, when 
conditions are streamline Re,/[k,/kg)* approaches 
the value for isotopes. The deviation produced 
by turbulence is apparent, but there is no 
attempt at a correlation. Similar conclusions 
apply to the results in Table 19. 


29-1 
270 


expt. 
T,°C 120 
AT 100 
ke 
kg x wt 70 
Re, 


Re. 
| 
K 


| 0-407 


A COMPARISON OF THE EXPERIMENTAL 
AND CALCULATED VaLtues OF K ann H 


The individual K and H values have been given 
in the preceding tables. These are now brought 
together, and their ratios considered in the 
following Tables 20 and 21; in this way a com- 
parison can readily be made. 


Discussion 


It will be apparent from Tables 20 and 21, that 
the values of H and K are affected very con- 
siderably by the withdrawal rate c. This problem 
does not occur in the design of columns for isotope 
separation, as the withdrawal rates are so small, 
that they do not interfere with the conditions in 
the annulus to the same extent. The considerable 
interference evidence with common gases is duc 
to the effect of convective turbulence Re, and 
turbulence produced by the mass flow resulting 
from continuous withdrawal, i.e. Re. The con- 
vective turbulence Re,, as indicated for the total 
reflux operation, may not be so adverse in effect. 
being counter-balanced by an effective reduction 
in annulus width leading to enrichment. The 
optimum value of Re,, and its relationship to the 
critical Re, of 25, will depend on each case, and 
has not been ascertained. 

Both for H, /CO, and H, /CO there is an increase 
in Hy,» and K,,,, with increase in o, depending 
on the AT values. However, no correlation has 
been attempted. 

It is apparent, that in the design of columns 
for separating common gases, pilot plant informa- 
tion is essential, The values of H and K must 
be available for the conditions anticipated for 
full-scale operation. Values of x would have to 


120 
270 
200 250 250 
5 5-90 | 18-50 0-43 
124 204 49 
885 BBO 160 21 
384 3-90 5-80 2-20 
| 1-83 | 3-66 | 6-96 5-91 | 6-55 18-52 0-455 
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Table 20. Column 2. H,/CO,; cy = 22-0% x = 0-25. 


H theor. K theor. 
H expt. Keheor. 


T = 50°C, AT = 60, [K/H),.%. = 505, 
0-0024 2-13 1-83 
0-0024 1-97 1-83 
0-0024 3-70 1-83 


T = 70°C, AT = 100, [K/H expt. = 368, 
0-0054 1-64 3-66 
00054 2-70 3-66 
0-0054 490 3-66 


T = 97:5°C, AT = 155, [H/K),.. = 276, Re, = 828 
0-0134 3-83 6-96 
0-0134 4-47 6-96 
0-0184 5-36 6-96 
0-0184 8-38 6-96 


T = 120°C, AT = 200, [K/H),..54. = 227 


0-0185 2-04 5-91 
0-0185 3-46 5-91 
0-0185 8-65 5-91 
0-0185 8-00 5-91 


T = 145°C, AT = 250, [K/H) 5.4. = 205 


0-0174 3-85 6-55 
0-0174 414 6-55 
0-0174 3-48 6°55 
0-0174 2-52 6-55 


be available for the mixtures being studied, and Hy.o Hex. for different temperatures, pressures. 
the influence of temperature, concentration, and and columns dimensions, coupled with a know- 
pressure on these values established. It is ledge of the flow conditions, such as Re, and Re. 
suggested that the ratios Kyyeor/Kexy., and could be used for design. 


Table 21. Column 2. H,/CO,.; cy = 15.0% CO. 
H theor. 
H expt. 
T = 70°C, AT = 100, [K/I1],,5:, = 360, Re, = 17-9 


0-00198 1:27 0-407 0-560 
0-00198 1-72 0-407 0-414 


H theor. K theor. K theor. 


T = 120°C, AT = 200, [K/H]),..4. = 247, Re, = 19-1 


0-00198 2-38 0-407 0-518 
0-00198 8-41 0-407 0-360 


T = 145°C, AT = 250, [K/H),..5. = 217, Re, = 28-5 
0-00198 1-98 0-407 0-636 
0-00198 2-90 0-407 0-425 
0-00198 7-58 0-407 0-163 


Re, = 47-4 
26-9 0-00113 0-57 3-21 210 
274 0-00121 0-61 3-00 150 
28-6 0-00065 0-326 5-61 
Re, = 66-5 
26-9 0-0083 1:225 2-98 630 
28-85 0-0020 0-728 5-08 250 
30-60 0-0011 | 0-408 8-98 70 
28°85 0-0085 7°25 670 
30-90 0-0030 8-33 390 
83-10 0-0025 10-10 160 
35-18 0-0016 16-30 20 
| 
82:40 0-0089 0-878 6-74 520 
85-50 0-00387 0-845 7-00 250 
88-00 0-0045 1-015 5-82 60 
29-90 0-0052 1-055 6-21 1220 
38-50 0-042 0-865 7:58 64) 
36-75 0-050 0-995 6-60 350 
39-80 0-0069 1-420 4-60 100 
theor. 
H 
20°45 0-00156 0-725 3380 
21-85 0-00115 0-084 140 
22-4 0-0021 1-09 54) 
24°85 0-00146 1-57 190 
22-05 000204 0-95 780 
26-00 0-00196 1-42 260 
28-85 0-00075 3-70 20 
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Comments on the operation of Columns | and 3 


According to Furry 
Haw -B kguw-B 


It is evident from these relationships that, for 
design, a compromise must be arrived at which 
will give a reasonable transport and enrichment. 

Column 1 has a large annulus, and conse- 
quently gives high transport and low separation. 
This was demonstrated. 

Column 3 has a small annulus, and the reverse 
effect is true. 

Column 2 was taken as the compromise. 

Column 3 is of greater interest than Column 1. 
as it gives high enrichment. From the results for 
the total reflux operation of Column 8, the gas 
mixture being H,/CO, (Table 2), it can be shown 
that Re = 1.86, k, = 0.00063, and k, = 0.00315. 
In this case, streamline conditions exist, and 
unlike the other columns, back-diffusion is greater 
than convective re-mixing. The asymmetrical 
term k, could not be calculated as the peripheral 
temperature variation of the inner tube could 
not be measured. The inner tube was centralized. 
The conditions could have been improved if « 
copper instead of a stainless steel inner tube had 
heen used. The value of k, must have been 
appreciable in this case, as an enrichment up to 
93-0°,, of CO, in H, at equilibrium was obtained, 
whereas according to the other conditions, a near 
100°,, separation should have been possible. 

To this effect of k, in column 8, would be 
added the effect of withdrawal rate, which would 
have a disastrous effect on the separation. Such 
columns could therefore be used only under 
carefully controlled conditions, and in « cascade 
system. 

CONCLUSIONS 

(1) The hydrodynamic equations, for discontin- 
uous operation of thermal diffusion columns, 
as deduced for isotopes, apply to common 
The adverse re-mixing effect of convective 
turbulence in the annulus is masked by what 
is probably a reduction in effective annulus 
width, resulting in increased separation which 
offsets loss of enrichment due to re-mixing. 


(iti) There is a lack of reliable information on the 
large effects of concentration and temperature 
on x values for various mixtures of common 
gases. This makes design of columns from 
first principles unreliable. 


From the results of continuous operation, it 
has been demonstrated that the values of H 
and K are dependent on withdrawal rates to 
such an extent that the theoretical approach 
to design as used for isotope separation is 
unreliable. 


Pilot-plant data are essential for mixtures of 
different gases to obtain column dimensions 
and optimum flow conditions for any desired 
separation. 

NOTATION 


L = heated length of column 
ideal separation length 
B ~ mean circumference of concentric tubes 
2u — annulus width 
‘9 = initial concentration of heavy componen 
¢, = final concentration of heavy component 
c = intermediate concentrations 
c, = equilibrium concentration of heavy compo- 
nent 
D, = coefficient of natural diffusion 
Dy = coefficient of thermal diffusion (truce) 
ky = thermal diffusion coefficient 
a = thermal diffusion constant 
7 = net differential transport 
H = transport coefficient 
Remixing coefficients 
ke = convective 
kq = back-diffusion 
kp = asymmetrical 
K = total 
N = number of ideal separation lengths in 
height L 
N = N/2-208 
Re, = Reynold's Number (convective) 
Re = Reynold’s Number of mass flow due to 
withdrawal 
T,°K = temperature of the cold wall 
T,°K = temperature of the hot wall 
T°K = arithmetic mean temperature across annulus 
p = gas density 
= gas viscosity 
o = withdrawal rate, ccs /min 
eo = withdrawal rate, grs/sec 
n =o /H = column factor 
Az, = column factor (vide Jensen [5] for defini- 
tion 
g = acceleration of gravity 
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Heat transfer to boiling water-methylethylketone mixtures 


A. S. Vos and S. J. D. van StraLen 
Laboratory of Physics and Meteorology, Agricultural University, Wageningen, The Netherlands 


Abstract—An apparatus has been designed for determination of heat transfer to boiling mixtures 
of water-methylethylketone under atmospheric pressure, using the central portion of a horizontal 
platinum heating wire at the same time as a resistance thermometer, according to the method 
of McApams 

A hysteresis effect observed in mixtures containing an excess of methylethylketone was avoided 
by keeping the annealed wire immersed in the liquid for a sufficiently long time at room tem- 
perature. The reference value of the resistance thermometer was aso measured with immersed 
wire at the boiling temperature of the liquid. 

The entire region of heat transfer up to the maximum of nucleate boiling was measured for 
some mixtures including the pure components (0, 42, 20, 52, 88-5, and 100% wt of methyl- 
ethylketone). With increasing concentration of methylethylketone a gradual! shift of the curves 
to lower heat transfer occurred, except for the 42% and 20% mixtures, where a noticeably 
high maximum heat flux of 2:5 and 2-0 times that of water was found. Since it was also observed 
with heating wires of different metals and alloys, this peculiar behaviour can be considered to 
be a characteristic property of the liquid mixtures themselves. 1956 


Résumé—Pour déterminer le coefficient de transmission thermique des mélanges d'eau et de 
méthyl-ethyl-cétone bouillant sous la pression atmosphérique, les auteurs ont concu un appareil 
utilisant la partie centrale horizontale d'un fil de platine servant d'élément chauffant ainsi 
que de thermométre a résistance selon la méthode de McApams. 

Ils ont évité I'hystéresis, observée pour les mélanges contenant un excés de methy! ethy! 
cétone en immergeant le fil recuit dans le liquide pendant un temps suffisant 4 la température 
ambiante. 

De méme, ils ont pris comme résistance de référence, celle du fil immergé dans le liquide 
bouillant. Pour des mélanges de composants purs (0, 4.2, 20, 52, 88.5 et 100% de methyl éthy! 
cétone), ils ont mesué les diverses valeurs du transfert jusqu'au maximum “ nucleate boiling." 

Pour une augmentation de la concentration de m.e.c., il se produit un déplacement gradue! 
des courbes vers un transfert plus petit, excepté pour les mélanges & 4,2 et 20% pour lesquels il 
y a un maximum notable du flux de chaleur, 2,5 et 2 fois celuide l'eau. Puisque ce 
particulier a été observé avec d'autres fils chauffants, en métaux et alliages divers, il s'agit bien 
d'une propriété caractéristique des mélanges liquides eux-mémes. 


METHOD AND APPARATUS could be used as a heating surface and at the 
The rate of heat flow g from a heated wire to a #me time as a resistance thermometer. 
boiling liquid is a function of 6 = t — T, the The heat flux © and the temperature differ- 
difference in temperature between the heating A 
ence 8 were directly calculated from measurements 


surface and the bulk liquid : 
=¢(8) a) of the potential drop E across the test section, 
— which had a length L, and of the current J, since 


The method described by McApams [5] has in 
principle been used to determine the heat transfer 
curves (1) for mixtures of water and methyl- 
ethylketone of different composition under atmos- 
pheric pressure. A horizontal platinum wire (of 
diameter D) was provided with thin potential 
leads, and thereby the central portion of the wire 


the specific resistance of platinum is known as a 
function of temperature [1]. For these measure- 
ments D.c. was used. One has : 

El 


q = 
(cal em”) = 0.07604 


(volts amperes cm™), and (2) 


_ 
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E 

= ¥(T) (3) 
Here A denotes the area of the test section, Ry 
the resistance of that section at the boiling 
temperature of the liquid in ohms, and ¥(T) a 
multiplication factor which for small or moderate 
temperature differences (i.e. in the ranges of con- 
vection and nucleate boiling) is found to be equa! 
to: 

1+ BT — CT 

B— 2CT 


Physically pure (99-99%) platinum wires, 
obtained from Drijfhout’s Edelmetaalbedrijven. 
Amsterdam, have been used, for which B is taken 
as 38-9788 x 10°. and C is taken as 5.88 x 107 
if T is expressed in degrees centigrade. The 
diameter of the heating wires was 0.020 and of 
the potential leads 0-005 cm. 

The chemically pure methylethylketone 
(M.E.K.) was obtained from Brocades-Stheeman 
& Pharmacia, Meppel, and distilled water was 
used. The important data pertaining to these 
mixtures were taken from literature [4], i.e. 


B.P. of methylethylketone : 79-56°C. 

Region of demixing : from 18-2 to 84.3% wt 
of M.E.K. 

B.P. in the region of demixing : 78-62°C. 

Composition of azeotrope : 88-6% wt of 
M.E.K. 

B.P. of azeotropic composition : 73-56°C. 


= (4) 


The test section of the heating wire was con- 
nected in series with a constant resistance of 
manganine. The potential drop across this 
resistance was directly proportional to the heating 
current. This voltage and that across the test 
section were recorded on a synchro-printing 
12-record fast speed Brown ElectroniK potentio- 
meter, with scale range from — 0-5 to 2-5 mV, 
after reduction to suitable values by application 
of tapped shunts. The reference value Ry was 
determined by passing a small current through 
the circuit while the wire was immersed in the 
boiling liquid. An anticorrosion circuit [5] was 
used by keeping the base of the boiling vessel at 
a tension 2 V lower than the negative end of the 
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Fic. 1. Diagram of boiling vessel. 
A = revolving arm; B = nickel coated brass base ; ( = 
total reflux condenser; D = draw bars; G = glass 
cylinder; M = calibrated thermometer; N = neoprene 
packing; L, and L, = potential leads; P = pertinax 
cover; S,. Sg, T, and T, = nickel coated brass bars ; 
W = heating wire. 


The boiling vessel (Fig. 1) consisted of a glass 
cylinder with ground endfaces. A nickel-coated 
brass cover ring and base were fastened by draw 
bars to the glass cylinder, and neoprene rings 
were used for packing. A pertinax cover was 
provided with four nickel-plated brass bars. One 
of the bars was fitted with a revolving arm, which 
made it possible to use wires of various length. 
The wire was stretched between this arm and the 
opposite bar, while the potential leads were 
attached to the other bars. A total reflux con- 
denser was mounted on the cover plate, and the 
boiling vessel was only half filled with the mixture. 
The liquid was then heated from below by a 
Bunsen flame until gentle boiling took place. 
This insured a uniform temperature throughout 
the liquid. 
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RESULTS WITH 
Piratinum Heatinc WIRES 


A curve giving heat flux 4 as a function of 


temperature difference © was determined by 
successively increasing the heating current, 
starting with a low value of J, as well as by 
decreasing the current starting with a high value 
of I. Two different curves were obtained, 
especially with mixtures in which the concen- 
tration of M.E.K. was high (Fig. 2- curve 2). 
This hysteresis effect was, however, removed 
(curve 8) when the wire, after annealing, was 
kept immersed in the liquid at room temperature 
for a sufficiently long period + (for which 18h 
was taken), before the observations were made. 
This procedure was adopted in al! further experi- 
ments, even in case of pure water (curve 1). 


Figs. 3 and 4 show heat flux 4, and co- 


efficient of heat transfer h = 4 respectively, 


as functions of temperature difference 6. With 
the same wire, which had a length L = 4-88 and 
a diameter D =0-0198cm, we _ successively 
determined the curves for : (a) water, (b) M.E.K. 
(8 x), (c) water, (d) the azeotropic mixture con- 
taining 88.5% wt of M.E.K. (3 x), and (e) a 
mixture containing 52% wt of M.E.K. The 
heating wire burned out after these determinations. 
and the mixtures (g) and (h), containing 20% 
and 42% wt of M.E.K. respectively, were 
investigated, using other heating wires. 

In the cases (a) to (d) inclusive, increasing 
current was first applied, and the current was 
reduced before the maximum heat transfer in 
nucleate boiling was reached. After these 
measurements, an entire curve was determined 
for increasing current only, until maximum 
nucleate heat transfer was obtained. Since the 


—= (BTU ft 2) 


Fic. 2. Effect of immersion of platinum wire on heat transfer for convection and nucleate boiling to water and 
to methylethylketone. 


The curves represent heat flux us u function of temperature difference between heating surface and liquid. 
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Heat flux for convection and nucleate boiling to water-methylethylketone mixtures as a function of tem- 
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Fig. 3. 
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results obtained with increasing and decreasing 
current were found to tally satisfactorily, a 
determination with increasing current only was 
made in the experiments (e), (g) and (h). 
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—= BFC) 
Fic. 4. Coefficient of heat transfer for convection and 
nucleate boiling to water-methylethylketone mixtures as 
a function of temperature difference between heating sur- 
face and liquid. 
The figures at the curves denote %, wt of M.E.K. 


In Figs. 2 and 8, all curves show a region in 
which convection seems to be rather inefficient 
(“ micro "-convection). Similar phenomena have 
been described by other investigators (8; 6; 7). At 
« temperature difference @ of about 10°C a 
stronger convection sets in. 

A gradual decrease of heat transfer occurs with 
increasing concentration of M.E.K., if this con- 
centration exceeds 52% wt. The nucleate boil- 
ing maxima of heat flux and of coefficient of heat 
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transfer in the 52% and 88.5% wt of M.E.K. do 
not exceed the corresponding value in water. On 
the contrary, however, in the 4.2% and 20% wt 
of M.E.K., the maxima exceed the value in 
water by not less than 150% and 100% in heat 
flux, and 70% and 85% in coefficient of heat 
transfer, respectively. It may be useful to 
observe that, for these mixtures, the size of 
vapour bubbles leaving the heating surface is 
definitely smaller than for the other compositions 
investigated. 


Resutts with HeatTinc WIRES or 
DIFFERENT METALS AND ALLOYS 


With the object of ascertaining that the peculiar 
behaviour of the heat transfer in the 4.2% and 
20% mixtures is really a property of the liquid 
mixture itself, and not a property of the special 
combination of heating surface and mixture, 
wires consisting of other metals have been used. 
These investigations have been carried out with 
a.c. and the maximum heat flux in the region of 
nucleate boiling only has been determined. 

The results are shown in Fig. 5. For these 
experiments, 20% and 28% wt of M.E.K. were 
used. The ratio of the maximum heat flux in 
these mixtures to that in water is also shown 
(Fig. 5- curves 1 to 10). An increase similar to 
that with platinum wires was observed in all 
cases, although the numerical value of the ratio 
varied. 

With a nichrome V (80 Ni; 20Cr) wire of a 
diameter D = 0-08 cm, which had been oxidized 
in air, the very high heat flux of 75 cal sec’ em™ 
in case of a 4-2% wt of M.E.K. was found. The 
maximum heat transfer to water was also con- 
siderably higher than with the other wires used, 
viz. 40 cal sec™* cm™ (Fig. 5 — curve 11). FARBER 
and Scoran [2] determined the heat transfer to 
boiling water with wires of a diameter of 0-1 cm, 
which were oxidized at red heat in boiling water. 
They found maxima of 47 cal sec"'cm™ for 
chromel A (80 Ni; 20 Cr), and an average of 
72 cal sec" cm™ for chromel C (57 Ni; 14-18 Cr; 
28-28 Fe) respectively, at atmospheric pressurc. 
Unfortunately chrome! C was not at our disposal, 
but it seems not unreasonable to expect that the 
42% mixture would show a still considerably 
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Heat transfer to boiling water-methylethylketone mixtures 


boiling heat flux as compared with the pure 
components must be expected, and at which 
concentration this should occur, will be con- 
sidered in a subsequent publication [8]. 


= 
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—= 10 (VA) max. (BTU hr-ttt-2) 
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NOTATION 


A = Area of surface of test section em? 


B,C = Constants in temperature equation 


2 
3 

4 5 

| : | | 

Fic. 5. Maximum heat flux for nucleate boiling of different h= 4; = Coefficient of heat transfer in test section 

wires to 42%, 20% and 23% wt of methylethylketone cal sec“! em= °C“! 


mixtures in comparison with maximum heat flux to water. I = Electrical current strength through test 
1. Pt (D =002cm); 2. Ag (0-02); 8. Cu (0-025); 4. Ni section A 


D = Diameter of test section em 
i E = Potential drop across test section v 


(0-02); 5. oxidized Ni (0-02); 6. Fe (0-025); 7. W (0-041); = Length of test section em 
8. constantan (0-02); 9. nichrome V (0-02); 10. oxidized = Heat transfer rate in test section cal sec" 
nichrome V (0-02); 11. oxidized nichrome V (0-08). == Sent Guu tm test caction 


cal 


higher maximum heat flux with this oxidized cesistance of test ot belling 
temperature of liquid Q 


alloy. = Temperature of test section 
From the results obtained with different heating T = Boiling point of liquid at atmospheric 

wires, it is concluded that high maximum heat pressure c 
transfer in nucleate boiling occurring in certain © = ‘— T = Temperature difference between test section 
mixtures of water and methylethylketone is a and belling Hquid, hes 

= Multiplication factor in temperature differ- 
characteristic property of the mixtures them- ence equation °C 
selves. The question whether for mixtures in + = Time of immersion of wire in liquid at room 
general an increase of the maximum nucleate temperature h 
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Book review 


Ullmanns Encyklopaedie der Technischen Chemie 
(8rd Edition), Vol. 5,854 pp. 237 illustrations. DM 108 
£9 9s. Od. (18 vols.). Urban & Schwarzenberg, Munich. 


Tue volumes of this monumental third edition of Ulimann's 
Encyclopaedia are being published in quick succession. 
So far the first volume of the general introduction to 
chemical engineering has appeared (the second volume 
is eagerly awaited), as well as three of the ten volumes 
which cover, in alphabetical arrangement, the whole 
field of applied chemistry including its theoretical founda- 
tions. 

The approach, in this volume which ranges from calcium 
carbide to di-isosyanate, is as comprehensive and lucid 
as in the previous volumes. References to pertinent 
literature and patents have been included up to 1953. 
The uniformly high standard of accuracy and modernity 
is remarkable, and the treatment seems to the reviewer 
free from any noticeable blemishes or errors. For those 
with sufficient knowledge of German this encyclopaedia 
is of great value. 


Announcement 


CHEMICAL ENGINEERING CONFERENCE 


An International Conference with the theme 
“ Chemical Engineering in the Coal Industry " is being 
sponsored this summer by the National Coal Board. It 
will be held from 26th to 29th June, 1956, at the new 
laboratories of the Board's Coal Research Establishment 
at Stoke Orchard, near Cheltenham, Gloucestershire. 


The Conference will comprise four technical sessions : 
“The Physical and Chemical Pre-treatment of Coal,” 
“ Fluid Bed Carbonisation,” “ Briquette Carbonisation " 
and “ Liquid By-Products." Delegates will receive and 
discuss seven papers given by English, French, German 
Dutch and American scientists and chemical engineers 
which will deal with the fundamentals and with the 
practical aspects of these subjects. English, French and 
German will be the working languages of the Conference. 
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